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Abstract 

The existence of a finite basis of algebraically independent one-loop integrals has underpinned 
important developments in the computation of one-loop amplitudes in field theories and gauge 
theories in particular. We give an explicit construction reducing integrals with massless propagators 
to a finite basis for planar integrals at two loops, both to all orders in the dimensional regulator e, 
and also when all integrals are truncated to 0(e). We show how to reorganize integration- by-parts 
equations to obtain elements of the first basis efficiently, and how to use Gram determinants to 
obtain additional linear relations reducing this all-orders basis to the second one. The techniques 
we present should apply to non-planar integrals, to integrals with massive propagators, and beyond 
two loops as well. 

PACS numbers: 11.15.-q, 11.15.Bt, 11.55.-m, 12.38.Bx 
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I. INTRODUCTION 



The computation of higher-order corrections to amplitudes in gauge theories is important 
to the experimental program at particle colliders. Recent years have witnessed dramatic 
advances in technologies for computing one-loop amplitudes, critical to the program of next- 
to-leading order calculations for collider physics. 

Important advances have also been made in computations of amplitudes beyond one 
loop. The computation of two-loop amplitudes relies in part on the ability to compute 
two-loop integrals, which has seen remarkable progress in the last decade. Several tech- 
nologies [1J] played a role in these advances, most notably_the Mellin-Barnes approach to 



computing integrals pioneered by Smirnov [2] and Tausk |3J, and later automated by Cza- 
kon [4]. Smirnov and Smirnov have recently introduced an alternative automated strategy 
for resolving singularities jsj. Anastasiou et al. have developed another method 6] of in- 
tegral evaluation combining sector decomposition [7] with contour deformation 8 . These 
technologies have played a key role in higher-loop calculations in the M = 4 supersymmetric 
gauge theory |9l. Il0|. 

The computation of amplitudes has also made use of techniques for reducing arbitrary 
tensor integrals to a basis set of scalar master integrals. In calculations performed to date, 
the reductions have relied on integration by parts (IBP) [ll|] to construct linear equations 
relating the various integrals, and on Gaussian elimination in the form of the Laporta al- 
gorithm [l^ to solve them. The solution determines a set of master integrals, and gives 
expressions for the remaining integrals in terms of them. This reduction approach has 
been automated in Anastasiou and Lazopoulos's air program 
gram 14|, and more recently, in Studerus's Reduze program [15|], as well as various private 



131 ]. in Smirnov's FlRE pro- 



computer codes. We should note that the existence of a method, such as the Mellin-Barnes 
approach, for evaluating loop integrals directly means that a reduction to master integrals 
is not, strictly speaking, necessary for a Feynman- diagram calculation. It greatly reduces 
the complexity and difficulty of such calculations, however. In order to use master integrals 
in such calculations, one needs the explicit forms of the reduction equations. 

Recent years have also witnessed the development and elaboration of a new set of tech- 



nologies, so-called on-shell methods 
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-|22j. for computing amplitudes. These rely only on 



knowledge gleaned from physical states. The unitarity method, one of the tools in this 
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approach, determines the rational coefficients of loop integrals in terms of products of tree 
amplitudes corresponding to cutting propagators in the loop amplitude. (These coefficients 
are rational in spinor variables.) 

It is possible to determine the set of loop integrals that contribute to a given process 
during the computation of their coefficients, and most of the higher-loop computations to 



date have proceeded in t 



generalized unitarity 



13, 



lis manner. In the most powerful form of the unitarity method, 



23H25], one cuts an amplitude into more than two pieces; indeed, 



in 'maximal unitarity' one cuts as many propagators as possible in a given contribution, 
thereby reducing any higher-loop amplitude to a product of basic tree amplitudes. The 
power of this technique is greatly enhanced by an a priori knowledge of a basis of integrals, 
as it then becomes possible to design the cuts in a general way. 

Knowledge of a basis, in contrast, is essential to developing an automated numerical 
implementation, which several groups are currently pursuing at one loop 26]. The required 
basis has been known for a long time at one loop. The four-dimensional one dates back 

n 

to the work of Melrose [27]. It is worth noting that the reduction equations themselves 
are not required when using generalized unitarity, because the method avoids the need 
for reductions of integrals with non-trivial numerators. We only need to know the set 
of algebraically independent master integrals. Baikov's work [28] suggests an interesting 
connection between integration by parts and maximal unitarity. 

Smirnov and Petukhov 29] have recently shown that the integral basis resulting from 
integration by parts is finite. In this paper, we give an explicit reduction to a finite set of 
integrals for planar integrals at two loops. There are two different kinds of bases we will 
consider. One requires algebraic independence to all orders in the dimensional regulator 
e (a "D-dimensional basis"), while the other requires algebraic independence for integrals 
truncated to O(e ) (a "regulated four-dimensional basis"). The latter contains fewer inte- 
grals, and is the relevant basis for the computation of amplitudes for numerical applications. 
We shall show how to limit the set of planar integrals that enter into a general two-loop 
computation, and discuss the reductions of some of these integrals. We leave the complete 
enumeration of basis integrals, as well as proofs of their algebraic independence, to future 
work. 

The approach we will pursue here makes use of a chosen subset of IBP equations, de- 
signed to avoid the introduction of unwanted integrals with doubled propagators, as well as 
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supplementary Gram-determinant equations to take advantage of additional reductions pos- 
sible when the loop integrals are performed in a truncated expansion about four dimensions 
rather than in arbitrary dimension (that is, to all orders in the e expansion). The approach 



we describe should also apply to non-p 



the Mellin-Barnes approach 



flfl 
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anar integrals, and beyond two loops as well. We use 



3l| to cross-check our equations, alongwith another 



technique for evaluating general higher-loop integrals, sector decomposition |7|, |31|. 

We will not discuss the analytic evaluation of the master integrals. Integrals involving 
a single dimensionless ratio of invariants may be expressed in terms of harmonic polyloga- 



rithms introduced by Vermaseren and Remiddi [32|, or alternatively in terms of the gener- 
alized polylogarithms of Goncharov {33)]; some integrals involving two dimensionless ratios 
can be expressed in terms of a two-dimensional generalization of harmonic polylogarithms 
introduced by Gehrmann and Remiddi |34j; for examples, see ref. |35j. The four-mass double 



box was computed by Ussyukina and Davydychev [36|, |37|. It is plausible that the complete 
set of two-loop basis integrals with massless internal lines can be expressed in terms of 
generalized polylogarithms, but this remains to be proven. 

In section HIl we review the basis of one-loop integrals with massless propagators in order 
to illustrate the two different bases, and to give a simple example of the use of Gram- 
determinant equations. In section 1111} we show how to reduce two-loop tensor and scalar 
integrals of sufficiently high multiplicity (again with massless propagators), thereby provid- 
ing a constructive demonstration of the existence of a finite basis. We also describe how to 
obtain a compact set of equations relating only integrals relevant to amplitudes, avoiding 
the introduction of integrals with doubled propagators. In section [V] we discuss the mass- 
less double box in detail. In section IVIj we apply these techniques to double-box integrals 
with different patterns of external masses. In section IVHt we apply the techniques to the 
pentabox integral. In section IVHH we give one example of the reduction of a six-point 
integral, the double pentagon. In section IIXI we present a heuristic explanation of some of 
our results using generalized unitarity. We summarize in a concluding section. 

II. REDUCTION OF ONE-LOOP INTEGRALS 

As a warm-up exercise, let us review integral bases at one loop along with their derivation. 
Throughout the paper, we will take the external momenta to be strictly four- dimensional. 
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(a) (b) (c) (d) 

FIG. 1: The basis of scalar integrals: (a) box, (b) triangle, (c) bubble, and (d) tadpole. Each corner 
can have one or more external momenta emerging from it. The tadpole integral (d) vanishes when 
all internal propagators are massless. 



They may be massless, or massive (representing, for example, sums of massless momenta 
in the original amplitude). In addition, we will take all vectors contracted with the loop 
momentum to be strictly four-dimensional as well. These vectors might be momenta or 
polarization vectors. All internal lines are taken to be massless. 

In an n-point one-loop amplitude in gauge theory, we start with integrals with up to n 
external legs, and up to n powers of the loop momentum in the numerator. (In a gravitational 
theory, we would start with up to 2n powers of the loop momentum. Up to questions 
regarding ultraviolet divergences, their treatment follows the same approach as the gauge- 
theory tensor integrals.) These powers are contracted with external momenta, external 
polarization vectors, or external currents. We shall denote the scalar integral by I n , 

I n (K l} ... } K n ) = J n [l] = -ij j—^ _ _ R ^ 2 — ^- — — — p . (2.1) 

In this equation, Kj...i = Kj + • • • + K\. We denote integrals with a function of i inserted in 
the numerator as follows, 

M' WJ - J (27r)D &{e _ _ . . . (£ _ KlHn _ 1}) 2 > M 

where the momentum arguments are left implicit. 

Let us begin with high-multiplicity integrals, with five or more external legs. Consider a 
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generic tensor integral 1 I n [£ ■ V\£ ■ v% • • • £ ■ v n ]. We shall make use of Gram determinants, 

g( Pi, "' ,Pi ) = det (2 Pi - qj ), (2.3) 
G( Pl ,---,Pi) = of* 1 ''"'* 1 ), (2.4) 



which have the useful property that they vanish if either the {pj} or the {qj} are linearly 
dependent. Using these objects, we can expand each of the four- dimensional vectors Vj in a 
basis of four chosen external momenta 61, 62, 63 > 64, 

1 



£(61,62,63,64) 



61, 62, 63, 64/ \6i, 6 2 , 63, 64/ V61, 6 2 , 63, 64/ 



,61,62,63,64 



(2.5) 



This leads us to consider integrals with numerator insertions where all of the Vj are equal 
to one of the 6j. At one loop, these factors are all reducible, because we can rewrite any dot 
product as a difference of denominators, for example, 

£ ■ 61 = \[{£ - Kf - {£ - K - h) 2 + (K + h) 2 - K 2 } . (2.6) 

The first two terms lead to integrals with fewer propagators and fewer powers of £ in the 
numerator, while the last two lead to integrals with fewer powers of £, 

I n [(t ■ v) n ) — ► I n -i[(£ ■ vy- 1 } © I n [(£ ■ v) n - 1 } . (2.7) 

Repeating this procedure for the daughter integrals (with a new basis element instead of 61 
where required) ultimately leads to integrals I n with n < 4 or with trivial numerators. 



n 



As is well known we can also reduce four- or fewer-point integrals with non- 

trivial numerators, by relying on Lorentz invariance to re-express them in terms of integrals 
where the non-trivial numerators involve only external momenta. (We could alternatively 
introduce additional basis vectors up to contributions from 0(e) numerators.) For the pur- 
poses of studying reductions, it therefore suffices to take v to be one of the external momenta 
even though the latter do not suffice to provide a basis. At higher loops, not all integrals can 
be reduced this way, because not all numerators can be written as differences of propagator 



1 While this numerator does not have free indices, they could be exhibited by differentiating with respect 
to the Vj, so in a slight abuse of language, we will refer to the integral as a tensor integral. 
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denominators as in eq. ( 12. 6p . Many integrals with irreducible numerators can nonetheless be 
simplified using IBP technology as implemented (for example) in the Anastasiou-Lazopoulos 



air code 



13| in Smirnov's Fire package [l4|, or in Studerus's Reduze package fl5|. leaving 



a smaller set of master integrals. 

We must next reduce the five- or higher-point integrals with trivial numerators ('scalar' 
integrals). While the reductions above hold independent of the dimensionality of the loop 
integration, the same is not true of all of the reductions we must consider here. We must 
distinguish between an integral basis to all orders in the dimensional regularization parame- 
ter e, and one which holds only through order O(e ). The latter may contain fewer integrals 
than the former, because it is possible for linear combinations of integrals to be nonzero 
but of O(e). At one loop, this is indeed what happens, with the scalar pentagon integral 
required in an all-orders basis while being reducible to O(e ) \^\ - 

Let first consider the reduction of six- or higher-point integrals, which can be done to all 
orders in e. Because the external momenta are four-dimensional, we have, 



G 



1,2,3,4 







(2.8) 



,5,1,2,3,4, 

where we have used the labels of the external momenta to represent the momenta themselves. 
Accordingly, 



G 



1,2,3,4 



,5,1,2,3,4 
If we now expand the Gram determinant, 







(n > 6) . 



(2.9) 



G 



£,1,2,3,4 
5,1,2,3,4 



G 



1,2,3,4 
5,2,3,4 



{i-K^G 



1,2,3,4 
5,^i 2 ,3,4 



\t-K 12 ) 2 G 



+{£-K 



1234 J 



G 



1,2,3,4 

5,1,^23,4 

1,2,3,4 

1,2,3,^45 



K 12 3YG^ 
— #12345) 2 G 



1,2,3,4 



5,1,2,^4 

1,2,3,4 
1,2,3,4 



-K(G[ /'J' 3 '/, ) +K 2 l2 G 



-#123 G 



+#12345 G 



5,^2,3,4 
1,2,3,4 

5,1,2,^4 

1,2,3,4 



1,2,3,4 

5,1,^23,4 



_ R 2 G , 1,2,3,4 
Kl2uU[ 1,2,3, K. 
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.1,2,3,4, 

we obtain an equation relating the n-point integral to six (n — l)-point integrals, 



(2.10) 



I n (K u . . . , K n ) = c 1 I n _ 1 (K nl ,K 2 ,...,K n _ 1 ) + c 2 In-i(K 12 ,K 3 



K n ) 



7 



+C 3 I n - 1 (K 1 , ^23, K 4 ,..., K n ) + C4l n -l(K U K 2 , ^34, #5, • • • , #„) (2.11) 
+C 5 I n -l(Kl, • • • , #45, • • • , Kn) + C 6 I n -i(K 1 , . . . , X 56 , . . . , K n ) 



where 



, 2 1,2,3,4 \ , - / 1,2,3,4 \ - ^ 1,2,3,4 



l ,5,iT 12 ,3,4y 12 U,l,^ 23 ,4y 123 V5, 1,2,^4 



2 , 1,2,3,4 \ 2 /1,2,3,4 
" Kl234 G 1 1, 2, 3, ^ 45 J + ^ 12345 G ll, 2, 3, 4 



1 ^/1,2,3,4 
Cl " ~o U,2,3,4 

C2 = -^^ 1 ' 2 ' 3 'A), (-) 



c V5,Xi2,3,4 

1 / 1,2,3,4 
c 3 = — G 



C \5, 1,^23,4/ ' 

1 r ( 1,2,3,4 
4 " c G U,l,2,K 34 

1,2,3,4 

C \1, 2, 3,-^45 

1 ^/1,2,3,4 
c 6 = — G 



c Vl,2,3,4, 

One can check numerically that the coefficient Co does not vanish for generic momenta, and 
hence the q are well-defined. 

In D dimensions, as mentioned above, pentagon integrals are needed as independent 
basis elements. When expanding about D = 4 — 2e dimensions, however, only the (9(e) 
terms are independent, so that the integral can be eliminated from the basis. We can 
derive the required equation by considering the Gram determinant G(£, 1, 2, 3, 4). The Gram 
determinant itself is of 0(e), because it can avoid vanishing only when the e components of 
£ appear in place of £. This leads us to consider the integral, 

h[G(£, 1,2, 3, 4)]. (2.13) 

One might worry that the Gram determinant can end up multiplying divergent terms in 
the integrand, yielding terms which are overall still of 0(e°) or even divergent. However, 
all divergences of the integral arise from regions where £ is soft, or collinear to one of the 
external legs. In these regions, the Gram determinant vanishes. Because the divergences are 
logarithmic at D — 4, any vanishing of the integrand suffices to eliminate the divergences. 
(At one loop, this in fact follows directly from the dependence of the integral only on the 
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e-dimensional components of t\ but that will not necessarily be manifestly true for similar 
integrals we shall consider in the two-loop case.) Furthermore, the integral is ultraviolet- 
finite by power counting. Accordingly the integral itself is also of 0(e), 

J 5 [G(£,l,2,3,4)] = 0(e). (2.14) 

We can use this to obtain a useful equation for the pentagon integral by expanding the Gram 
determinant, and re-expressing dot products of the loop momenta in terms of differences of 
denominators, 

G(£, 1, 2, 3, 4) = d + d x £ 2 + d 2 (£ - K x f + d 3 (£ - K 12 ) 2 + d 4 (£ - K l23 ) 2 + d 5 (£ - K 12U ) 2 
where 

2\2nf K l2, 3,4\ 2 2 fK 12 ,3,A\ 2 2 ^fl,K 23 ,4 



+ 2^^ 3lG (^ 2 3 3 4 ) - 2^ G ( 1 #3 4 ) H- 2^ 3 ^ 3lG ( 1 - 2 2 ^) 



(A'J 231 ) 2 G'' L2 



1,2,3 



4 = 2 G ( 1 < 2 ' 3 '^ -kIg( k ^ 4 ) + k1 2 gC' k ^) -KMC**' 

VI, 2,3,47 1 V 2,3,4 J 12 V 2,3,4 7 123 V 2,3,4 



d 2 = kIg( k ^ a ) -k 2 2 g( k ^) + k 2 123 g( k ^/) -k 2 12U g( Ki ^ A 
1 \k 12 ,?>a) 12 u,K 2 3,4y 123 u,2,K 34 y 1234 v 1,2,3 

^3 = -k 2 g( k ^ a ) + k 2 2 g(^ k ^ a )- 

1 \l,K 23 ,Aj 12 U,^ 23 ,4y 123 \1,2,K 3 J 1234 V 1,2,3 

di = kIg( k ^ a )-k1 2 g^ k ^ 

1 VI, 2, AW 12 U,2,AW 123 VI, 2, K 3 J 1234 V 1,2,3 



W A 12 ,3,4\ , - /1,/fas, 4\ . /l,2,lf M \ , r ^ 2 ^M,3 



5 1 V 1,2,3 / 12 V 1,2,3 / 123 V 1,2 , 3 / 1234 VI, 2, 3 

The integrals of the terms on the last two lines of eq. (I2.15P will vanish, as they correspond 
to box integrals with e(£, . . .) in the numerator. The integral of the do term is simply a 
pentagon, and the integrals of the di,...,5 terms are box integrals. These reductions and 
relations yield the well-known basis shown in fig. [TJ 
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Our aim is to extend these considerations to two-loop integrals. We will also introduce 
a new technique for reducing integrals with irreducible numerators to the set of master 
integrals, based on rewriting the system of IBP equations. We delineate the finite universal 
basis, while leaving a complete and detailed enumeration of it to future work. 

III. REDUCTION OF PLANAR TWO-LOOP INTEGRALS 
A. The Integrals 




(c) 

FIG. 2: The three basic types of two-loop planar integrals, labeled by the number of legs attached 
to each leg of the vacuum diagram: (a) P ni ,ri2i (b) -Pni,n 2 » ( c ) Pn*n 2 - 

We turn now to our main object of study, the planar two-loop integrals. We can organize 
the different integral skeletons we obtain, representing only the propagators, into five classes. 
Three classes, depicted in fig. [2} arise from attaching external legs to the non-trivial two-loop 
vacuum diagram shown in fig. [3J We obtain planar integrals by attaching external legs to one 
or two of the internal lines, and possibly to its vertices. Were we to attach external legs to 
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FIG. 3: The non-trivial two-loop vacuum diagram. 





(a) 



(b) 



FIG. 4: Two-loop integrals which are products of one-loop integrals, labeled by the number of legs 
attached to each leg of the vacuum diagram: (a) I ni>m , (b) I* ± m . 

the third internal line as well (here, the middle line), we would obtain non-planar integrals. 
This gives rise to three of the five types of two-loop planar integrals; the remaining types 
are simply products of one-loop integrals. We label the integrals according to the number of 
external legs attached to each of the vacuum diagram's internal lines, denoting the absence 
of lines attached to vertices by stars. The three types of integrals are, 



ni,n 2 



d D t x d D e 2 1 

(2^(2^ Wi - ^i) 2 • ■ • - K^Wi + ?2 + K ni+n2+2 y 



ni,n 2 



£ 2 2 (£ 2 -K ni+n2+1 ) 2 ---(£ 2 -K {ni 



(ni+2)-(ni+n 2 + l), 



d°£ X d D £ 2 1 



(2tt)^ (2tt)^ - K x y ■■■(£,- K^yfr + £ 2 f 

1 

£ 2 {£2 — K ni+n2 + i) 2 ■ ■ ■ (£ 2 — K( ni+ 2)...( ni +ra2+l), 

d D £ x d D £ 2 1 



(2n) D (2n) D £\{£ x - K x ) 2 • • • {h - K x ..^)\l x + £ 2 f 

1 

V 

£\[£ 2 - K ni+n2 ) 2 ■■■{£ 2 - K {ni+1) ... {ni+n2) ) 2 
along with products of two one-loop integrals shown in fig. HJ 



—i 



., / d D £ x d D £ 2 1 



(2vr)^ (2tt)^ l\{i x - K x f ■■■(£,- K x ... ni f 



(3.1) 
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v 

4(^2 — -^ni+n 2 +l) 2 • • • (4 — -^(ni+2)---(n.i+n 2 +l)) 2 

d D l x d D £ 2 1 



' " ' (2n) D (2n) D f 1 (£ 1 -K 1 y---(t 1 -K 1 ... ni y 

1 



(3.2) 



£ 2 (£ 2 - K ni+n2 ) 2 • • • (4 - i^(„ 1+ i)...(„ 1+ „ 2 )) 2 

so that P ni> „ 2 is an (ni+n,2+2)-point integral, P* n2 and J nijn2 are (n!+n 2 +l)-point integrals, 
and -P^* n2 and i£ ina are (ni + n2)-point integrals. Without loss of generality, we may take 
Tlx > n 2 . 

In our discussion below, we will focus on the P ni ,n 2 integrals. Similar arguments typically 
apply to the P£'*n 2 - We may also observe that, 

^ Z> ni,n 2 (-^1 ? • • • > -^ni+n 2 +l) Pn\,ri2\K\i • • • i -^ni+na+1) 0) 5 

-fni.^i-^ 1 ' ■ ■ ■ ' -K-m+ni) = P^^i^i, ■ ■ ■ , K ni ,0, K ni+ i, . . . , K ni+n2 ) . (3.3) 

so that the values of P* and P** are known in terms of P. Nonetheless, their different 
branch cut structures strongly suggest that the former are algebraically independent of the 
latter. In explicit examples in sections FvTlVHIl we examine P** integrals. 

These integrals will arise in the leading-color contributions to two-loop QCD amplitudes, 
including amplitudes for production of electroweak bosons or other particles coupled to 
quarks. Just as in the one-loop case, all internal lines will be massless, but the external 
legs of the integrals can correspond to sums of external momenta, and hence can be either 
massless or massive. Each of the vertices can come with a power of the corresponding loop 
momentum, and each of the three-point internal vertices in P^*<** can also come with a 
power of either i\ or £ 2 , so that we should consider tensor integrals with up to (rii + 2) 
powers of l\ along with n 2 powers of £ 2 , or alternatively up to (ni + 1, n 2 + 1) or [n\^n 2 + 2) 
powers of the two loop momenta (£i,£ 2 ). 



B. Reduction of High-Multiplicity Integrals with Non- Trivial Numerators 

We begin our discussion of integral reduction at two loops by considering tensor integrals 
with rii > 4, P nitn2 [£ ■ v\ £ ■ v 2 ■ ■ ■ £ ■ v n ], where each £ can be either £\ or £ 2 . We can use 
the expansion of eq. (12. 5ft . with the external momenta bx, . . . , 64 chosen amongst the first n\ 
momenta. This leads us to consider integrals with numerators containing factors of £\ ■ Kj, 
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where 1 < j < n\. As in the one-loop case, these numerators are reducible, 

li ■ Ki = \[t!i ~ Ki-u-i)) 2 - {h - K lmj f + Kl.. 3 - Kl. u ^] . (3.4) 

The first two terms lead to integrals with smaller indices (P ni _i in2 , P^-i „ 2 , or one of P niiTn _i 
and P* j in the case n\ = n 2 ) and simpler tensors, while the last two lead to integrals 
with simpler tensors. Repeating this procedure, including application to leads to tensor 
integrals P ni ,n 2 with m < 4 and n 2 < 4, or integrals with general (ni,n 2 ) but with trivial 
numerators. 



C. Reduction of High-Multiplicity Integrals with Trivial Numerators 

Once we have eliminated high-multiplicity tensor integrals, or equivalently those with 
non-trivial polynomials in £ 12 in the numerator, we must consider integrals with trivial 
numerators but arbitrary number of external legs. To reduce P ni ,n 2 integrals with n\ > 5, 
we can make use of the same Gram determinant as in the one-loop case, 

flf'^fUo. (3.5) 



5,1,2,3,4 



We will obtain the reduction, 



-fnijTiaV-^lj • • • > -^ni+ri2+2) 

Cl-Pni-l,n 2 (-^2) • • • ) K(n 1 +n 2 +2)l) + C 2-Pni -l,n 2 (-^12 > -^3; • • • > K ni+n2+2 ) (3-6) 
+C3P ni -l,n 2 (Kl, K23, K±, • • • , K ni+n2+2 ) + C 4 P ni _i in2 (i^i, K 2 , -^34, K5, . . . , i^„ 1+n2+2 ) 
+C5P ni -l,n 2 (-^l) • • • > K45, • • • > ^m+712+2) + CeP ni -l,n 2 (-^1 j • • • ) -^56; • • • ; -^ni+n 2 +2) , 

where the coefficients q are given in eq. (I2.12p . the same ones as in the one- loop 
reduction. This reduction involves only propagators dependent solely on l\. (In 
the case ri\ = n 2 , P„ 1 _i in2 ({i^ i }™ 1 ; {i^}^^ 2 " 1 " 1 ) is given by the flipped integral 
Pn 2 ,ni-i({-^}nJ+n 2 +ii {-^i}ni)0 This reduces the set of scalar integrals to P ni ,n 2 with 
n>2 < ni < 4. This means we have a finite (if large) set of integrals in terms of which 
we can express any planar two-loop integral, and hence any planar two-loop amplitude. 

This reduction generalizes both to non-planar and to higher-loop integrals. The details 
are in some cases more intricate, but at two loops, we can reduce all integrals with more 
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than eleven propagators. We can see this using a variety of Gram determinant equations 
similar to eq. ( I3.5|) . 

There is another, more general, way of looking at this question. Let us label the momenta 
in the two-loop vacuum diagram of fig. E]by £ 1; £ 2 , and £3. They are not independent, because 
^1 + ^2 + ^3 = 0, but can be treated symmetrically. There can at most be eleven independent 
invariants ti involving the loop momenta, namely the three squares of the loop momenta, 

£l 4 £l; (3.7) 

and eight invariants built out of loop momenta, of the form 

lj ■ h , (3.8) 

where four hi are selected out of the external momenta. Because the external momenta 
are strictly four-dimensional, we can express all remaining ones, and hence all invariants 
involving them, in terms of these four. We can choose the invariants to be all manifestly 
reducible. If we have more than eight external lines attached to the lines in the vacuum 
diagram (that is, excluding legs attached directly to the vertices in it), then there are more 
than eleven propagators with denominators di, and accordingly we can write down non- 
trivial equations, 

= di-J2cjtj (3-9) 

j 

for denominators beyond the eleventh. Inserting this equation into the integrand allows us 
to reduce the integrals with more than eleven propagators to simpler integrals, at arbitrary 
D or equivalently to all orders in e. (This assumes that the coefficient of the original integral 
is non-zero, which the earlier discussion demonstrates for the two- loop case.) 

D. Integration-by-Parts Without Doubled Propagators 

The reductions discussed in the previous subsections show that all required integrals for 
a planar n-point two-loop amplitude can be written in terms of the integral with a 
trivial numerator; P^ n2< ^ integrals with trivial numerators or numerators dependent only 
on £ 2 ; and -P n i<4,n 2 <ni integrals with trivial numerators or numerators dependent on either or 
both of the loop momenta. Of course, some numerators can still be written as a difference 
of denominators, as in eq. ( 13. 4p . The corresponding integrals can then be reduced. The 
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remaining irreducible integrals, for n% + n 2 > 5, are those with irreducible numerators, 
which cannot be written in such a way. For example, in P^, £ 2 ■ K\ would be an irreducible 
numerator. (Some integrals with n\ + n 2 < 4 require a more specialized analysis, just as at 
one loop, and for the most part we shall not consider them in the present article.) 

To reduce these integrals, where possible, we will employ the integration-by-parts (IBP) 
technique first introduced long ago by Chetyrkin and Tkachov and refined into a general- 
purpose algorithm by Laporta |l2| . 

The IBP technique as outlined by Laporta, and as implemented in air jl^ j. Fire jl^ . 
and Reduze I 



P 



a ti a i 2 o v 

J (2n) D (27i) D dt.D(LJoAKA) ' 1 ' ; 



relies on writing down all possible equations from introducing a differenti- 
ation inside the integrand, 

,.l r d D i x d D £ 2 d 

(2n) D (2n) D dFj D(£ u £ 2 , {K t }) 

where D(£i, £ 2 , {Ki}) is the denominator found in eq. ( 13. ip . The simplest choices for are 
the set of external momenta, along with the two loop momenta. The use of dimensional 
regularization ensures that there is no boundary term in this equation. 

With these choices, however, the resulting equations involve not only the integrals of 
interest (as well as simpler planar integrals), but also integrals with doubled propagators. 
These arise from derivatives hitting the denominators of the integrals. Moreover, such 
integrals can have worse infrared singularities, so that their use results in the appearance 
of additional inverse powers of e in coefficients. This in turn would require them to be 
known (either analytically or numerically) to higher orders in e. These integrals do not 
arise directly in the computation of gauge-theory amplitudes, and usually we do not wish 
to introduce them at the stage of solving these equations. In all cases studied to date, it 
has been possible to eliminate such integrals (at a cost of retaining some integrals with non- 
trivial numerators), and it seems plausible that this holds true generally. Their elimination 
requires the consideration of very large systems of equations, performed using the Laporta 
algorithm. In addition to the considerable computational complexity of these systems, which 
has made it difficult to proceed in explicit examples beyond four-point integrals, it is also 
far from clear how to characterize these systems in the general case. 

For this reason, we seek to simplify the equations we obtain by eliminating unwanted 
integrals, those with doubled propagators, from the very start. We will do so by making 
special choices of the vectors in eq. (I3.10p . For example, we could choose vectors whose 
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dot product with the numerator resulting from differentiating any propagator vanishes, 

v(£-K) = 0. (3.11) 



As these expressions are the coefficients of the doubled propagators, this vanishing will 
ensure that doubled propagators are absent. We can construct such vectors using Gram 
determinants; defining 

q( V, K 62, 63, h \ = _^_ G f w > & i> & 2, h, h \ , 3 12 . 

\di, d 2 , d 3 , d±, d$) dwa \di, d 2 , d 3 , d 4 , d 5 J ' 



^ 6, 7,8- 



vectors t> of the form 

v» = G( ti > t 

\£ 2 , £ u 1,2,3,4, 

will have the desired property with respect to propagators depending only on l 2 . For exam- 
ple, the IBP equation, 

P^[d 2 -v\ = (3.14) 

(where o> 3 - = d/d£j), will be free of doubled propagators. Because £^ 2 are D-dimensional 
vectors, Gram determinants like that in eq. (13 . 13[) give the most general solution to eq. (13. lip . 

However, this solution is not general enough for our purposes. The problem is that while 
the constraint (13. lip is sufficient, it is not necessary. It is in general too strong a constraint, 
and in practice we would miss equations if we insisted on it. The weaker constraint, that 

v ■ (£ - K) oc {£ - Kf , (3.15) 

suffices to remove the doubled propagator as well. This constraint can be expressed as 
the requirement that there be no remainder upon synthetic division of v • (£ — K) by the 
propagator denominator (£ — K) 2 . We must impose this constraint for every propagator. 
For an integral with rid propagators, that is with a denominator in the form, 

W- 1 = II d s = n(*ii'i + ^£2 - K 3 ) 2 , (3.16) 

3=1 3=1 

we must impose the rid equations 

Rem 11 *f 311 J2 2 ^ = 0, (3.17) 

(anh + a i2 £ 2 - KA 2 



or equivalently, 



Rem [^i + ^]-(^ 1 + ^ 2 -^-) = 
(ajtit + a j2 £ 2 - KjY 
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where Rem denotes the remainder on synthetic division (using either £j as a variable). (In 
these equations, the <ij will be ±1 or 0.) These equations are for vectors v\$ built out of 
the loop momenta, external momenta, and dot products thereof. We will discuss how to 
find the general solution to these equations in the next section. (By convention, we order 
the denominators as follows: first those depending only on the first loop momentum; then 
those depending only on the second loop momentum; and finally, those depending on both 
loop momenta.) 

Before trying to solve the equations, let us try to characterize the solutions better. These 
have several general properties that will be helpful in finding and using these vectors. For 
example, if we have a pair of vectors, {vf\ v^} that satisfy eq. (I3.17p . then any multiple 
of the pair is also a solution. In particular, multiplying by any Lorentz invariant involving 
either of the two loop momenta gives us a solution. 

Not all these additional solutions are useful, however. We can divide these Lorentz invari- 
ants into two types: the reducible ones, expressible as a linear combination of propagator 
denominators and external invariants, and irreducible ones whose dependence on the loop 
momenta cannot be expressed using propagator denominators. While multiplying by an 
invariant of the former type does yield a solution to the constraints (I3.17p . it is not a useful 
solution, because it does not yield an independent equation for the integrals of irreducible 
numerators. To see this, let us write out the resulting IBP equations. The original equation 
is, 

• (vxW) + d 2 ■ (v 2 W)} = . (3.19) 
Multiplying the vectors by a factor / gives a sum of two terms, 

I[f(9i ■ (viW) + d 2 ■ (v 2 W))} + I[W(v x ■ d x f + v 2 ■ d 2 f)} = . (3.20) 

The first term contains reduced integrals (that is, with fewer propagators) and terms pro- 
portional to the original equation. What about the second term? If / is reducible, its 
derivative can be written as a linear combination of the derivatives of the denominators 
in W. Because of eq. A3. 18j) . the sum in parenthesis is strictly reducible, that is reducible 
without adding any terms proportional to external invariants. The second term in eq. (13.201) 
thus contributes only reduced integrals. Accordingly, only irreducible factors / can give rise 
to new IBP equations. 
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Indeed, the solution itself won't be useful if both of the vectors v\ and v 2 are reducible; 
we will therefore restrict attention to solutions in which at least one is irreducible, that is, 
at least one term in one of the pair t> 12 is irreducible. Note that not all independent vectors 
will lead to independent IBP equations; but because the independence of the IBP equations 
can depend on whether the dimensional regulator e is taken to zero or not, we leave that 
assessment to a later stage. 

It is possible to find even weaker constraints that remove double propagators, by adding 
a 'total derivative', that is a function which integrates to zero, to the right-hand side of 
eq. f l3.15p . We will not consider such right-hand sides. In the examples we consider below, 
they are not necessary, though we know of no general proof. 

Intcgration-by-parts equations can be supplemented with Lorentz-invariance equa- 



tions 



34|, using operators built out of derivatives with respect to external momenta. In 



general, these are not independent of the complete tower of IBP equations [42] ; because we 
will be able to generate the complete tower of IBP equations, we do not need to consider 
Lorentz-invariance equations. 

As mentioned in the introduction, for generalized unitarity, we only need to know the 
basis integrals, that is the set of integrals left independent by the full set of IBP equations. 
In order to find this set, we could of course solve for the IBP-generating vectors analytically, 
and then construct the set of IBP equations analytically as well. However, it suffices to 
solve for these vectors for a randomly-chosen ('generic') numerical configuration of external 
momenta. For higher-point integrals, or integrals with many massive external legs, this can 
greatly reduce the complexity of the calculation, and in particular the memory required to 
solve for the IBP-generating vectors. 



E. Additional Identities to 0(e) 

The integration-by-parts identities give relations between different integrals that are valid 
for arbitrary dimension D, or equivalently, to all orders in the dimensional regulator e. 
However, in practical calculations at a given order in perturbation theory we are interested 
in computing terms only through O(e ), and we are quite willing to drop terms of 0(e) or 
higher. Additional relations between integrals, even if they are only valid through O(e ), 
are for these practical purposes just as good as relations that hold to all orders in e. The 
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reduction of the one-loop pentagon integral, as we reviewed in section [TTJ, is exactly this kind 
of relation. In general, we can write down several forms of integrands leading to integrals of 
0(e), built of Gram determinants or products thereof, 



G(£ u b u b 2 , 6 3 , b 4 )G(£ 2 , b[, b 2 , b 3 , 6' 4 ) ; G 
r (£i,i2,bi,b 2 ,b 3 \ r (£i,£2,bi,b 2 ,b 3 ,b±\ 

where the momenta attached to the first loop (through which l\ flows) are contained either 
within the set {b±, b 2 , b 3 } or within the set {b'{, b 2 , b 3 }, and the momenta attached to the sec- 
ond loop are contained within the other of the two sets. These Gram determinants all vanish 
when either loop momentum approaches a potential (on-shell) collinear or soft configuration, 
thereby removing the corresponding divergences from the integral, and rendering it finite. 
In addition, the Gram determinants vanish when both loop momenta are four- dimensional, 
so that the integrals are of 0(e). We can also write down differences of expressions yielding 
finite integrals which will again vanish when both loop momenta are four- dimensional, so 
that the resulting integrals are again of 0(e), 

£i,h,b 2 ,b 3 \ f£ 2 ,b[,b 2 ,b 3 \ fl u b u b 2 , b 3 \^f£ 2 , b[, b> 2 , b' 3 



\£xM,b 5 ibJ Wb'^b'J \£ 2 ,b',ib' 5 ,b'J WhMibe) 1 ' ) 

where the legs attached to the first loop are all represented amongst the 6j, and the legs 
attached to the second loop, amongst the b\. (For P* and P** integrals, k\ must also be 
amongst the b[, and k n amongst the b{, for P**, k ni+ \ must be amongst the b i} and k ni 
amongst the &■.) 

Not all of these determinants will necessarily lead to useful equations reducing the basis. 
We can also consider integrals with numerators containing a product of one of these Gram 
determinants and other irreducible factors, so long as the integrals are ultraviolet-finite 
(which can be determined by power-counting). As is true for the IBP generating vectors, 
we can also generate additional identities for a randomly-chosen configuration of external 
momenta; this will be sufficient to identify the integrals that are independent through O(e ). 
Typically, we will first solve all D- dimensional IBP equations, and use the solutions of 
those equations (in analytical or numerical form) to reduce the integrals obtained from 
inserting Gram determinants into the numerator; this will provide additional identities to 
O(e ) between the independent master integrals. 
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We can write down additional Gram determinants beyond those given in eq. (13.21 j) . 



■1,^2, h,b 2 ,b 3 \ ^f^kMMM 



u^b'^b'j' u 2 ,b[,b' 2 ,b' 3 ,b'J ' 1 • j 

where all momenta attached to the second loop (with loop momentum £ 2 ) are represented 
amongst {b\,b 2 ,b 3 } in the first case, and similarly for the momenta attached to the first 
loop in the second case. However, these determinants give rise to integrands which are in 
fact total derivatives, and hence the corresponding integrals vanish identically. To see this, 
consider the following vector, 

G [ e n u u u u > { 6 - M ) 



&2, &3> &4/ 

where all momenta attached to loop in which £ 2 flows are in {61,2,3 }• The vector's dot product 
with the derivative of any propagator with respect to £ 2 will vanish so that only the derivative 
of the Gram determinant itself can enter any equation; but that derivative is proportional 
to the first determinant in eq. (13. 23ft . Accordingly, we do not need to consider the forms in 
eq. (I3.23P if we have already solved the IBP equations. If we include additional irreducible 
prefactors, we will again either obtain an expression proportional to the determinants in 
eq. (I3.23P or to linear combinations of them and the last determinant in eq. (13. 21 p . Equations 
similar to those considered in this section were obtained for six-point integrals by Cachazo, 

n 

Spradlin, and Volovich [10J using leading singularities. 



IV. IBP-GENERATING VECTORS 



In order to find the general form of vectors leading to IBP equations free of doubled 
propagators, we must find the general solution to the set of equations (13 . 1 8D . We begin by 
rewriting them in a somewhat more convenient form, 

[a^vt + a j2 v 2 ] • {<Tj\£.i + a j2 £ 2 - Kj) + Uj {<J jX l\ + a j2 £ 2 - Kj) 2 = , (4.1) 

where Uj is a polynomial in the various independent Lorentz invariants of the loop and 
external momenta. Because the different propagator denominators are independent (the 
integrals for which this is not true we have already treated in sects. IIII Bl and IIII CI) , this 
equation must hold for each of the rid propagators independently. 
Let us also write a general form for the , 

«T = c? 1 ^ + c^t 2 + £ Ci ( V , (4.2) 

beB 
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where the sum runs over a set of n B — up to four — basis vectors for the external momenta. 
(There would be four basis vectors for integrals with five or more external legs, and n — 1 
vectors for integrals with fewer.) Each of the coefficients cf^ is again a polynomial in the 
various independent Lorentz invariants. 

We consider as independent variables only invariants that are independent with respect 
to the loop momenta. That is, i\ and i\ are independent, as are each of these with respect to 
£i-k\, and a given invariant of the external momenta, say k\-k 2 - However, different invariants 
of external momenta are not independent, which is to say their ratio should be treated as a 
constant parameter. Let us pick the one independent invariant to be s\ 2 = {k\ + k 2 ) 2 , and 
define the ratios, 

Si 



Xij 



svz 

S{.. .j 



Xi-j —■ (4-3) 

2 

mj 

Hi =—, 

Sl2 

in order to express the remaining invariants in terms of si 2 . (For certain integrals with fewer 
than four external legs, we should pick a different invariant.) We will term these quantities 
parameters. Each of the coefficients cf would have an expression in terms of the invariants, 

V = {£l,h- £ 2 , l\ , {h ■ b} beB , {£ 2 ■ b} beB , s 12 } . (4.4) 

We treat these invariants as the basic symbols or variables out of which we build the solutions. 
For example, coefficients of engineering dimension two could be expressed as follows, 

c? = c%s 12 + £ ■ b + E • b + c%el + cg^ • £ 2 + c%£j . (4.5) 

beB beB 

The coefficients c[ P j of each term are rational functions of the parameters Xi—j an d Hi- (In 
order to distinguish the different dimensions to which we will refer below, we refer to the 
engineering or energy dimension of £i and b as such, dropping the 'engineering' qualifier only 
when context makes it unnecessary.) 

Our discussion generalizes in a straightforward way both to higher loops and to non- 
planar integrals. At higher loops, we will have a vector vi for each loop; the expansion (14. 2 h 
will have a sum over all loop momenta; and the set of variables V in eq. (14.41) will include 
all squares of loop momenta, all dot products of loop momenta with each other, and all dot 
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products of the loop momenta with the basis vectors in B. In general, some dot products 
of loop momenta with other loop momenta will be irreducible, but this does not change 
the derivation of the equations. For non-planar integrals at two loops, we will have more 
than a single equation involving two different vectors. (This will anyway be true at higher 
loops.) The general approach to solving the equations we now outline will also carry over, 
though the specific procedure that solves the equations most efficiently in these more general 
settings remains to be investigated. Internal masses will introduce additional parameters in 
eq. (14. 3 p while leaving the basic invariants (I4.4p and the general structure unchanged. 

Without loss of generality, we consider only solutions V\ t 2 of homogeneous engineering 
dimension. We could in principle proceed by writing down a general form for coefficients of 
a given engineering dimension, starting with dimension zero, and proceeding by increments 
of two. Plugging in these general forms into eq. (14. 2p . and requiring that the coefficient 
of each monomial in the basic variables (14. 4p vanish independently, we would obtain the 
solutions of the given dimension. 

This method of solution works quite well for finding solutions of low dimension, but be- 
comes very memory-intensive for higher dimensions. Furthermore, it does not allow us to 
determine when we have found the complete independent basis set of solutions, namely the 
set of solutions t>i,2 to eq. f 14 . 2 [) in terms of which all others can written as linear combina- 
tions, with coefficients that are polynomials in the basic variables ( 14 . 4 j) and rational in the 
parameters (14.31) . 

As a reminder, we are interested only in solutions for which at least one coefficient in 
either v\ or t>2 is irreducible with respect to the set of denominators. Let us denote the 
operation of removing terms proportional to a propagator denominators (that is, reducing 
by the set of propagator denominators) by the operator Irred. It leaves behind only the 
irreducible part of an expression. (This operation is most naturally implemented using a 
Grobner basis for the propagator denominators, but as these are all linear in the basic 
variables (I4.4p . the use of such a basis is not essential.) We defer a precise definition to later 
in this section. 

We can assemble the set of equations eq. ( 14. ip into a single matrix equation. To do, first 
assemble the various coefficients using the relabeling, 
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and the definitions, 



= 2, 

pUA) = J + 2, j e {l,...,n B }, 

p(£ 1; 2) = n B + 3, (4.6) 

p(£ 2 ,2) = n B + 4, 

p(j,2) = n B + j +4,j e {l,...,n B }; 

c p{q ,i) = cf\ g G{4,4}U5,ze{l,2}, 
C 2 n s +4+j = Uj , j G {1, . . . , n d } . (4.7) 

Treating the coefficients c as a row vector 2 , define 

-E2n s +4+j,i = -dj • (4.8) 

S is a (2ns + 4 + n^) x nd-dimensional matrix; the number of rows we will label n r . Each 
column corresponds to eq. (14. ip for a different propagator. We then have the following 
matrix equation, 

cE = 0. (4.9) 

Mathematicians call each solution to this equation a syzygy of E. 

In intermediate stages, we may need to solve not only homogeneous equations such as 
this, but also inhomogeneous equations, 

cE = f, (4.10) 

where the row vector / is independent of c though it may depend on other parameters. 

Both of these equations are linear polynomial diophantine equations. In the adiatretoflu- 
ous language of mathematicians, the former is an equation for the syzygies of the ideal 
submodule of Q{{x%-ji l JL i\)W] nd generated by the rows of E, 

Syz(( ei ,...,e n ». (4.11) 



2 Row vectors provide a more natural interpretation, as this choice also leads to treating the derivatives of 
the equation with respect to this vector's entries as row vectors, which in turn leads to a more natural 
implementation in a symbolic algebra language such as Mathematica. 
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More precisely, we seek a linearly-independent basis for the irreducible elements of the syzygy 
module, 

Syz((e 1; . . . , en))/ Syz Red ((e!, . . . , e n )) , (4.12) 

where the R e d subscript denotes the reducible subspace with respect to Irred. In this lan- 
guage, it is basically a textbook problem, though there are aspects which require a bit more 
work than a textbook solution. 

The solution relies on the use of Grobner bases 43) . The reader may find an explanation of 



the varied uses of Grobner bases, as well as algorithms for their construction, and the required 



background material, in several textbooks 
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of that Adams and Loustaunau 



451 ]. Of these, we shall primarily make use 



44| . Sturmfels gave a brief overview of Grobner bases 46] , 



and Lin et al. ]47] also provide a nice introduction to Grobner bases of modules from a 
physicist's point of view. Grobner bases have been studied for use in integral reductions by 



Smirnov and Smirnov 



481 ]. and have been used by Smirnov in Fire [14], as well as in other 



studies of integral reductions [49 1. 

Our review here will mention only the minimal material needed for the description of the 
solution. Grobner bases, amongst other uses, provide a certain generalization of linearly- 
independent bases to nonlinear multivariate polynomials. The basic setting is that of poly- 
nomials in a set of symbols. In our context, the symbols are those in the set V (14.41) ; the 
coefficients are arbitrary rational functions of the Xi-j an d forming the field Q{{Xi -ji 
in mathematicians' language. We will need to consider vectors or tuples of polynomials as 
well as polynomials themselves. 

The basic machinery requires us to choose an ordering of the terms built out of the basic 
symbols, as well as of tuples of terms. There are various ways of doing this, of which a 
lexicographic ordering is conceptually the simplest. While the choice of ordering will not 
change the space of solutions we find, the efficiency of the (standard) algorithms we employ 
will depend greatly on this choice. We choose the so-called degree-reverse lexicographic 
order (DRL or grevlex) for the basic symbols, and a term-over-position (ToP) ordering for 
tuples of polynomials. A generic term in a polynomial built out of the symbols Xi has the 
form, 

cx v = cx p ^ •••<"> (4.13) 
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while a generic n-tuple of polynomials is a sum of terms of the form, 



cv = cx v e v = cxi 1 ■ ■ ■ x^ n e v , (4-14) 

when the unit basis tuples are the set {e,,}. The DRL ordering for polynomial terms starts 
with a basic ordering of the symbols x n , . . . , x\ and orders x p before x q (denoted x p -< x q ) 
iff, 

X Pi < X 9» or X/ Pi ~ X/ 9* an d ^ ne rightmost non-zero entry in p — qr is positive 

i i i i 

(4.15) 

The ToP ordering orders tuples containing a lone monomial as follows, 

v -< u x v -< x u or x v = x u and e v -< e u . (4.16) 

(The basis vectors ej are ordered by their first non-zero component.) The leading monomial 
of a polynomial (or tuple) P, denoted lm(P), is the monomial v which is last in the ordering, 
v y v' for all monomials v' in P (stripped of any coefficient c). 

With an ordering chosen, we can define a polynomial reduction algorithm, essentially a 
repeated synthetic division with respect to a basis set B of polynomials (or tuples), yielding 
a set of coefficients c, and a remainder r, 

p-^r, (4.17) 

where 

p=Y,Cbb + r. (4.18) 

beB 

The coefficients are again polynomials in the basic symbols, and the remainder is a poly- 
nomial (or tuple of polynomials if p is a tuple). At each stage of the synthetic division, a 
polynomial is divisible by a selected divisor iff its leading monomial is divisible by the divi- 
sor's leading monomial. For the purposes of synthetic division, we can treat n^- dimensional 
tuples of polynomials by taking their dot product with an n^-tuple of dummy or "tag" 
variables (t\, . . . ,t nd ), and then performing ordinary synthetic division with the set of vari- 
ables extended to include the tag variables t«. The ToP/DRL ordering is then given by a 
DRL ordering, with the tag variables ordered before the other variables. One must ensure 
that other algorithms used maintain the linearity in these tag variables. At the end of a 
calculation, the tuples can be recovered by differentiating with respect to them. 
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In general, the reduction coefficients c& in eq. (I4.18P are not universally defined; they 
depend not only on the ordering chosen for monomials, but also on the order in which 
the polynomials are taken during synthetic division. Only if the set B of polynomials is 
a Grobner basis will the reduction coefficients be independent of the order in which the 
polynomials are taken. 

Let us write out the method of solution we have used, reverting to physicists' language, 
and postponing until the next section an explicit example, that of the massless double box. 
We use the Buchberger algorithm to compute the required Grobner bases, though more 



sophisticated algorithms 50j are available and would be worth investigating. Using the 
algorithm, we compute the Grobner basis of the set of rows of E, treated as n^-tuples. We 
assemble the elements of the basis, again n^-tuples, into a matrix G. The number of rows is 
determined by the number of tuples n g in the Grobner basis, which may be smaller, equal 
to, or larger than the number of original vectors n r (which is equal to 2ne + 71^+4 in the case 
of E). In addition to the basis itself, we will need the cofactor matrix C, which expresses 
the basis elements in terms of the original vectors, 

G = CE. (4.19) 

It may be computed as a by-product of Buchberger's algorithm (or other algorithms) for 
computing the Grobner basis. Because G is a Grobner basis, we may also express each 
original vector as a linear combination of the basis vectors; this defines another matrix Q, 

E = QG. (4.20) 

In order to find the syzygies of E, we must first write down those of G. The syzygies of 
the latter — that is, n„-tuples s such that sG = — can be constructed as described in 



the textbooks [44|, |45[. The construction starts with the S-polynomial of two rows or basis 
n d -tuples g itj , 

e („ „\- MMgj>M n km(]m(gj),1m(9j)) „ a01 x 
bp{9i,9j) ~ lm( ^ 9i lm(f?j) Bj, (4.21) 

where 1cm denotes the least common multiple, with the added definition that lcm(ej, ej) = 
if i ^ j. (The factors in front of are then pure polynomial terms, with the basis vectors 
e canceling out.) The S'-polynomial also plays a central role in the Buchberger algorithm 
itself. By construction, this S'-polynomial can be completely reduced over the Grobner basis, 

n g 

S p (g i ,g j ) = £,hk9k-> ( 4 - 22 ) 



fc=i 
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Each syzygy or linear relation between the n g Grobner basis elements can be represented by 
an n g -tuple, with basis elements ct\, . . . , cr ng , 

<r k G = g k . (4.23) 

We can define a basic syzygy or linear relation, 



Si. 



lcm(lm(flfi),lm(gj)) lcm(lm(gi), lm(g 3 -)) ^ {j 
hnTO MgJ) ^ " g h * ° k ■ (4 ' 24) 

The complete set of syzygies of G is then generated by the set 

{% 1 1 < i < j < n g } . (4.25) 

These syzygies, treated as row vectors, can again be assembled into a matrix S. With S in 
hand, the rows of SC are syzygies of E. In addition, because we must have 

E = QG = QCE , (4.26) 

then (/ — QC)E = 0, where / is the identity matrix. The rows of / — QC are thus also 
syzygies of E. We extend SC to add these rows. In our application, these turn out to be 
relevant only in some variants of the solution algorithms described below. 

Not all syzygies, that is rows of 5*, are linearly independent. Indeed, there are typically 
dozens or even hundreds of syzygies. The set of syzygies can be reduced in a variety of ways, 
of which the two principal ones we use are polynomial reduction and numerically-assisted 
row reduction. Furthermore, as discussed earlier, we are interested only in rows of SC which 
are independent after reduction with respect to the propagator denominators, and this allows 
for additional reductions in number. For this purpose, we use the Irred operator described 
earlier. We define it as the polynomial reduction (element-by-element) with respect to a 
Grobner basis Gd of the propagator denominators, built over the symbols in V (here using 
a plain lexicographic ordering), 

p— ^Vlrredp. (4.27) 
We apply polynomial reduction to the rows of S, treated as n^-tuples, using a special ordering 



of the underlying variables, called the Schreyer ordering 44 , 5lJ . It is defined by 



x p ei -< y Q ej <^=^> \m{x p gi) -< \m(y q gj) or \m(x p gi) = lm(y q gj) and j < i . (4.28) 
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It is useful because the syzygy generators (14.251) form a Grobner basis with respect to this 
ordering. When reducing syzygies, we start with those of lowest engineering dimension, 
removing those which reduce to zero (and hence are linear combinations of other syzygies), 
and proceed incrementally in the engineering dimension. 

To find the set of fully-independent solutions modulo reducibility, we can proceed as 
follows to convert it to a linear algebra problem. We form the irreducible part of the 
solutions, 

and convert them into "tagged" polynomials using "tag" variables as described earlier. We 
now construct a vector space, in which each coordinate corresponds to a different monomial, 
and where each monomial that may appear in any of the tagged polynomials, or in any 
product of an irreducible polynomial times a tagged polynomial, is assigned a coordinate. 
Each tagged polynomial P (that is, each solution s) may then be mapped to a vector 
Vec(s), whose entries are rational functions of the x%j an d the \x\. Independence can then be 
determined by linear algebra (e.g. row reduction). We can check it numerically, by evaluating 
the solution for a given numerical choice of external momenta. For a given solution s, we also 
need to generate the vectors corresponding to multiples of s by a factor x built out of the 
variables in V. We can do this either by mapping the multiple, Vec(a;s), or by multiplying 
Vec(s) by the appropriate matrix. After removing linearly-dependent solutions, we usually 
end up with only a handful of independent syzygies, which we assemble into a matrix S. 
The general solution to eq. ( I4.9j) can then be written as follows, 

c=(pi ■■■Pn s )S, (4.30) 

where n s is the number of independent solutions, that is the rows of S; and where each pi 
is an arbitrary polynomial in the variables in V. Because the Grobner basis is finite, and of 
finite engineering dimension, the complete set of solutions is generated by a finite and finite- 
dimensional set, and any effectively reducible integral will be reduced by an IBP equation 
built using one of the basis elements (with a possible irreducible prefactor determined by 
the dimension of the numerator in the integral). 

In order to solve inhomogeneous equations such as eq. ( 14.1 Oft , which will arise in some of 
the variants of the solution algorithm presented below, we must first reduce the right-hand 
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side / over the Grobner basis of E, 



f = q f G + r f . 



(4.31) 



If rf is non- vanishing, the equation has no solution. In our context, rj will typically have one 
or more free coefficients pi (arbitrary polynomials), and we will need to impose additional 
constraints on them to ensure that rj becomes reducible over G. (It is not necessary to 
make it vanish strictly.) We can do so by solving the homogeneous equation, 



where g is an n 5 -tuple of dummy coefficients. The solution will express the free coefficients 
Pi in Tf in terms of a more constrained set p\. The expression in terms of p\ will now be 
reducible over G. Once we have solved this subsidiary equation (or if r/ vanishes to begin 
with), a particular solution to the inhomogeneous equation is given by 



because then q/CE = qjG = f. In our solutions of inhomogeneous equations, we will not 
be interested in the general solution, but it can be obtained by adding an arbitrary solution 
to the corresponding homogeneous equation with / set to zero (obtained following the steps 
discussed above). 

A. A Simple Algorithm 

The steps described in the previous section can be summarized in a simple algorithm 
(Algorithm I), which starts as input with a matrix E as in the form (14.91) . 

1. Compute the Grobner basis G and the cofactor matrix C for the set of n^-tuples given 
by the rows of E; 

2. Build the set of syzygies S of G using eq. (14.24)) ; 

3. Reduce the set of syzygies by synthetic division with respect to previous retained 
syzygies, discarding those with no remainder. It is best to proceed incrementally in 
the syzygies' engineering dimension; 



r f -gG = Q 



(4.32) 



c = q f C, 



(4.33) 
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4. Construct the matrix Q which expresses E in terms of G; 

5. Construct the set of solutions, SC along with / — QC; 

6. Reduce to a set of independent solutions with respect to reduction by the set of 
propagator denominators. 

This algorithm works nicely and quickly for simple cases, such as the massless double box 
discussed in more detail in the next section. However, it suffers from very memory-intensive 
(and slow) intermediate stages for more complicated cases such as the four-mass double 
box or the pentabox, when the number of x an d A* parameters grows. There is room for 
improvement, because a great deal of unnecessary information (pertaining to fully-reducible 
solutions to the equations) is computed in intermediate stages. 



B. An Improved Algorithm 

For these reasons, we use a somewhat more involved procedure. The greater complexity 
of the procedure is balanced by simpler execution at each stage. The basic idea is to split 
up the solution into several stages. We can split the matrix E and the desired coefficients c 
into reducible and irreducible parts, 

E irred = Irred E , 

girrcd = Irredg? (434) 

^nred j^j ^rirred 



i 



grcd = ~ _ girrcd _ 

At the first stage, we solve the homogeneous set of irreducible equations, 

girrcd^irrcd = q _ ^35) 

The full equation can then be rewritten as follows, 

? Gd E = _c irred £ red , (4.36) 

which is an inhomogeneous equation for c red in terms of the (now-known) irreducible poly- 
nomials c irrcd . It turns out to be better to solve these equations in two stages, first only the 
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rows of E arising from propagators involving only l\ or £2, but not both; and adding in the 
full set of equations in a second stage. 

In order to solve the inhomogeneous equations, we write out an auxiliary set of equations 
which impose reducibility on each of the coefficients c red , 



The sum could also be taken over the Grobner basis elements used to define the Irred oper- 
ator. There are advantages and disadvantages to simply adding these equations as auxiliary 
equations, with additional unknowns c a j, as opposed to substituting these expansions into 
the inhomogeneous equations (I4.36p . and both the memory usage and the time required 
for solving can depend sensitively on this choice. In the examples we have considered, it 
appears better not to substitute in the first of the two stages of solving the inhomogeneous 
equations, and better to substitute in the second stage. 

At the first stage, we select the columns in eq. ( I4.36P corresponding to propagators 
involving either t\ or £ 2 , but not both; in the planar case, this means all but the last 
column. Call the solutions of these equations c^. At the second stage, we split the c a , 



and solve eq. (I4.36P for <5c red , with all but the last column of the right-hand side replaced by 
zero, and with an additional reducibility constraint of the form ( I4.37P imposed on 5c red . 

The strategy for solving the inhomogeneous equations ( I4.10p can be summarized in the 
following algorithm (Algorithm II), 

1. Compute the Grobner basis G of the rows of E (it is sufficient to compute a partial 
Grobner basis limited to the maximal engineering dimension found in the right-hand 
side /), along with the cofactor matrix C; 

2. Reduce the right-hand side / over G, yield coefficients qj and a remainder r/; 

3. If the remainder ry of this reduction does not vanish identically, solve the equation 
rj — gG = with dummy coefficients g using Algorithm I (in practice, it is best to 
impose an engineering dimension limit in the intermediate steps of Algorithm I, and 
increment it until the solution converges to a stable one); 



3=1 



(4.37) 




;rcd 



(4.38) 
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4. If a dimension-limited Grobner basis was computed in step 1, and r/ was not identi- 
cally zero at step 3, extend the Grobner basis to the new (typically larger) maximal 
engineering dimension in the constrained form of /. Rather than starting from scratch, 
one can start from the original Grobner basis G, in which case the full cofactor matrix 
will be given by the product of the new and old matrices, Cf — C 2 Ci. 

5. The solution is then qfC. 

The improved strategy for solving the original equation ( 14. 9 p can then be summarized in 
the following algorithm (Algorithm III), 

1. Compute solutions to eq. (I4.35P using Algorithm I; 

2. Solve the rows in the inhomogeneous equation (14 .361) corresponding to propagators 
containing a lone loop momentum, for the reducible terms c red , along with the con- 
straint equations ( I4.37P expressing reducibility, using Algorithm II; 

3. Write each coefficient c red as a sum of this solution, and another coefficient <5c red , 
as in eq. (I4.38p . Solve the inhomogeneous equation corresponding to the propagator 
containing both loop momenta, along with constraint equations of step 2 with their 
right-hand sides set to zero, using Algorithm II (here it is better to substitute the 
reducibility constraint equations (14.371) back into the inhomogeneous equations); 

4. Reduce to a set of independent solutions with respect to reduction by the set of 
propagator denominators. 



V. THE MASSLESS DOUBLE BOX 



1 4 

FIG. 5: The double box P 2 ** 2 . 
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As an example of how to apply the ideas presented in the previous section, let us examine 
the planar double box, P 2 * 2 , with all external legs taken to be massless. The integral is shown 
in fig. The D- dimensional reductions were worked out several years ago for the integral 
with all external masses vanishing, using air [131 ] . The same reductions have been worked 
out (though not reported explicitly) for configurations with one external mass 34|. The 
three- and four-mass cases have not been worked out previously. We discuss the massive 
cases in the next section. 

As described in the previous section, we start by looking for vectors v± 2 that give rise 
to IBP equations free of doubled propagators. For the double box, we choose £1,2,4 as basis 
momenta, so the general form (14. 2 p becomes, 



c i c l i c i c 2 i c i ft l + c i ^2 i c i ft 4 ) 



.(1) 



(2), 



(5.1) 



where each of the coefficients is itself a function of Lorentz invariants in V, which here 
we can take to be the following set, 
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?i , £i ■ £ 2 , i\ , t\ ■ h , ii ■ k 2 , t\ ■ h , £ 2 ■ h , £ 2 • h , £2 • h , s 12 } . 



(5.2) 



The other dot products of the loop momenta can be expressed in terms of these via mo- 
mentum conservation, and the other independent invariant is treated as a multiple of s± 2 , 
s i4 — Xu s i2- There are two irreducible numerators, l\ • fc 4 and £ 2 -k\. 
The matrix E of eq. (14. 8[) then takes the form, 



E = 8 



if 

h-li 
h-£i 
h A -£ x 











£\/4 








V 




-k v £ 2 +£ v £ 2 

k 2 -£i-s 12 /2 
h-£i ~ X14S12/2 
















- k 2 -£ x + 
k 3 -£ 2 + k A -£ 2 + £ x 



kv 
k 2 - 
k A - 



\ - s 12 /2 
\ - sia/2 
9 i + W2 























£v£2 

l 2 

k v £ 2 

—k\-£ 2 - k r £ 2 - k 4 

h-£ 2 





£1/4 + S12/4 - k\-£\j2 - k 2 -£ 1 /2 











£1/4 






(5 
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£i + £ 1 -£ 2 

£ r £ 2 + £ 2 2 

fci-^i + ifei-4 

k 2 -£ x -k x -£ 2 -k 3 -£ 2 -k±-£ 2 

A: 4 -£! + fc 4 -£ 2 

-U-£ x + £ v £ 2 k r £ 1 + k 2 -£ 1 + £ 1 -£ 2 £\ + £v£ 2 

-k 4 -£ 2 +£j -k r £ 2 -k 4 -£ 2 + £ 2 £i-£ 2 + £ 2 2 

k v £ 2 - X14S12/2 k v £ 2 + s 12 /2 k 1 -£ 1 + k v £ 2 

(l+Xu)si 2 /2-k r £ 2 -k 3 -£ 2 -h-£ 2 s 12 /2-k v £ 2 -k 3 -£ 2 -k r £ 2 k 2 - £ x -k x - £ 2 -k v £ 2 -k 4 - £ 2 

h-£ 2 k 4 -£ 2 -su/2 k 4 -£i + k 4 -£ 2 









-k 4 -£ 2 /2 + £ 2 /A 

£ 2 /4+s 12 /4-k r £ 2 /2-k 4 -£ 2 /2 

£?/4 + £ r 4/2 + ^/4 



while the vector c is 



~ r _ ( r (ti) (ta) (1) (2) (4) tt a ) (1) (2) (4) v /, 

C — C x C x Cj C x C 2 C 2 C 2 C 2 C 2 K.1...7J . ^D.4J 

The leftmost three columns correspond to equations for the left-hand loop; the following 
three columns, equations for the right-hand loop; and the last column, an equation for the 
common propagator. 

The propagator denominators are, 

£\ ■ {h - hf ; (£1 - K l2 f ■ £\ ■ {£ 2 - h) 2 ; (£ 2 - K u ) 2 ; {£ x + £ 2 f . (5.5) 

Reducibility of an expression may be determined by reducing over a Grobner basis of these 
denominators, which using a plain lexical ordering is, 

£i-£ 2 , £\ , h-£ 2 , k 3 -£ 2 - si 2 /2 , £\ , k 2 -£ x - sn/2 , k v £ x ■ (5.6) 

The Irred operator gives the remainder after this reduction; for example, 

Irred(ai£^ + a 2 £\-£ 2 + a 3 £\ + a 4 £\-k\ + a 5 £i-k 2 + a e £i-k 4 + a 7 £ 2 -ki + ag£ 2 -k 3 
+a 9 £ 2 -k 4 + ai si 2 /2) 

= (a 5 + a 8 + a l0 )s 12 /2 + a 6 k A -£ l + a 7 kx-£ 2 . (5.7) 

When using algorithm I from the previous section to solve eq. (14. 9p . we first obtain a 
Grobner basis for the seven-tuples making up the rows of E. There are 50 tuples in this 
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basis, which in turn give rise to 167 syzygies, which we can represent as 17-tuples. Polynomial 
reduction and removal of completely- reducible syzygies (with respect to the basis in eq. (15.61) ) 
leaves us with 119 syzygies of the Grobner basis. These in turn give rise to 101 solutions of 
eq. ( I4.9p . of which three are independent. The matrix I — QC gives no additional solutions 
(which can be understood on dimensional grounds here). 

There is one solution whose coefficients are of engineering dimension two, 

v x , x = -2(* 4 -4 + l\)K ~ + (2*r4 - + (4* r 4 + 2/s 2 .4 + 2* 4 -4 - s X2 )e x , 
v x , 2 = 2{£\ - k 4 -£ 2 )k1 + t*k% + (2k v £ 2 + £%)K + (2*3-4 - 2*r4 - s 12 )t 2 ; (5.8) 
and two solutions with coefficients of engineering dimension four, 
v 2 . x = (-4* 2 -4*4-4 - 4*3-4^ + 4*4-4^ - 4*4-4^ - 4£ 2 x £ x -£ 2 - 2l\£\ - 2 Xx d\s x <i)K 

+(4* r 4* 4 -4 - 2k v e x el - 2k 2 -i x i\ - Ak 3 -£ 2 £l - Ak 4 -i 2 i\ - 444-4 - n\t\ 

+2£\s X2 - 2 X u£ 2 x s X2 )k^ + {-±k x -£ x £\ - 4* 2 -4^ + 2(t\) 2 + 2£\s X2 )k'i 
+(4* r 4*2-4 + 4(* 2 -4) 2 + 8* r 4*3-4 + Sk 2 -£ x k 3 -£ 2 + 8* 2 -4*4-4 + 8k x -£ x k 4 -£ 2 
+8* 2 -4*4-4 - 4* 4 -4^i + 8*i-44-4 + 8*2-44-4 + 4* r 4^ + 4* 2 -44 

-4*i-4si2 - 6* 2 -4*i2 - 4* 3 -4*i2 - 2* 4 -4*i2 - 4* 4 -4*i2 + l\s X2 + 2xi 4 4*i2 
-44-4*12 -2l\s X2 + 2s\ 2 )t x , (5.9) 

v 2 . 2 = (4*1-4*4-4 + 4*3-4*4-4 + 4*4-4*4-4 - 4* 4 -44-4 + 4* 3 -4^ - 4* 4 -4^ 

+64-4^2 + 4(^) 2 " 24-4*12 - 2xi44-4si2 - 2£\s X2 )k^ + (4*1-4*4-4 

-4*4-44-4 + 2* r 4^2 + 2*2-4^2 + 4*3-4^2 + 64-4^2 + 4(^) 2 - 2X144-4*12 

-24s 12 )k% + (4* r 4^2 + 4*2-4^2 + 44-4^2 + 2(^) 2 - 2£\s X2 )k1 

+ (-4*3-4*4-4 - 4(* 4 -4) 2 + 2*3-4^2 + 2*4-4^2 + 2*r4*12 - 2X14*3-4*12 

-2xi4*4-4*i2 + ^2*12 + 2xi4^ 2 *i2Ki + (4* r 4*r4 + 4*1-4*2-4 + 4*1-4*3-4 

+4*2-4*3-4 + 8(* 3 -4) 2 + 8*1-4*4-4 + 8*3-4*4-4 + 8*3-44-4 - 2* r 4^ 2 

-2*2-4^2 + 8*3-4^2 + 4*4-4^2 - 2*1-4*12 - 2Xl4*r4*12 - 2*2-4*12 
-2X14*2-4*12 - 8*3-4*12 + 2*4-4*12 ~ 4* 4 "4*12 ~ 64"4*12 ~ 4Xl44"4*12 

-2£ 2 2 s X2 + 2s 2 X2 )£ 2 l ; 

and 

„ 3il = - *iM _ m? _ <®I _ 4 ^. fe _ m± _ 2k2 . eA _ m 

Xl4 Xl4 Xl4 Xl4 
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OK .f f 2 

+X1AS12 ~ 2£ v £ 2 s l2 - 2 X u£r£ 2 s 12 )k» + (4^-^4-4 - 2k r £ 2 il - 1 

+2jfe 3 -^? + 2k,-i 2 i\ " 2(^) 2 - ^ - 4^ r £ 2 - + 2kvi.il ~ 2l\l\ - 

+Xi4^si2 - 2 X idi-hsi 2 )k^ + (-4fci^ 2 A; 2 .4 + 2k v i 2 i{ + 2fc 3 -^ + 2£; 4 -^ 2 

Xl4 Xl4 

-4A: 2 -£ 1 fc 3 -^ 2 - 4/51-^4^1 - 4A; 3 -£ 2 A;4-£i - kk 2 -£ x k^£ 2 - kk±-i x k±-i 2 + 4^-^ 

Xl4 Xl4 Xl4 

, 4fc 2 -£A-£ 2 2fc r ^ „„ 2 2hriA 9 , „ , 
H 1 h 4/c 2 -£i£ 2 H 2xuki-iis 12 - Ak v £ 2 s 12 

Xl4 Xl4 Xl4 

2&r4si2 



2xi4^2-^iSi2 + 2/c 3 -£ 2 si 2 + 2xi4^3^2Si 2 + 2k 4 -£ 2 s 12 + 2xuh-£ 2 s 12 
Xu 

-2£\S 12 - £ -bl- 2 l±±^l + 4*14*1 - 4*12 - ^12 - ^ + X14S? 2 )^ + (2*1-^12 
Xl4 Xl4 Xl4 

+2xi 4 A:i^iSi2 + 2xi 4 A; 2 -£iSi 2 - 2A; 4 -£ lSl2 - 2^s 12 - 2x14^12)^ , (5.10) 

( Ai DO* 2k,-t 2 t\ Ah-i^-i, n2 2k x -£ 2 £ 2 2 
^3;2 = (-4fc 4 -£ 2 q 8k4-£ 2 h-h |-6fc r £ 2 £ 2 H h2fc 3 -£ 2 £ 2 

Xl4 Xl4 Xl4 

, „ fl2 2A; 4 -£ 2 £ 2 2£ 2 £ 2 4£ r £ 2 £ 2 2 (i 2 ) 2 

-2ki-i 2 i\ =— + U\i\ + — ^ + 8i x -i 2 i 2 2 + 2 + 2(£ 2 ) 2 + v 2; 



Xl4 ~ Xl4 " Xl4 ~ Xl4 

-2£ 2 s 12 - 3xi4^ 2 si 2 )^ + (4k r £ 2 k 4 -£ 2 + 4h-£ 2 £ 2 2 + 2 - 2k r £ 2 £\ - 2k 4 -£ 2 £ 2 2 

Xl4 

£ 2 £ 2 9ii -tJ 2 ( f 2 \ 2 

+2 £ 2 £ 2 + bh + U v i 2 i\ + + 2(£ 2 ) 2 + m. - 3 Xl4 f 2 s 12 )^ + (4(A:r£ 2 ) 2 

Xl4 Xl4 Xl4 
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+4fc r ^ 3 .£ 2 + 4^.^4-4 + 4A; r £ 2 £ 2 + + 8k 1 -e 2 e 1 -e 2 + 4fcr4£l " £2 

Xl4 Xl4 

+ 2k v £ 2 £\ + _ 2jfe 3 .^l - 2ife 4 .^l + 2£ 2 £ 2 + 111 2 + 4£ r £ 2 £ 2 + flli 2 ^ 2 

Xl4 Xl4 Xl4 

+2 (£ 2 ) 2 + (®! + 2^ Sl2 )^ + (-4(A; r 4) 2 + 4fcl ^ 2A:3 " £2 - 4(A; 3 -£ 2 ) 2 + 4^-^4-4 

Xl4 Xl4 

4A;i -4^4-4 „, „ , „ , , 2 2A; r £ 2 £ 2 „ n n2 2k v £ 2 £\ 
+ 2 -4A; 3 -£ 2 A; 4 -f 2 -4fc r £ 2 £ 2 ^^+4fc 3 -£ 2 £ 2 + 



fi , , ±k x -l 2 l x -l 2 , ou e e e , 4fc 3 -^i^ 2 2fc r f 2 f 2 . , 

-8k 1 -£ 2 £ 1 -£ 2 h 8k r i 2 i 1 -i 2 H + 4k 3 -£ 2 £ 2 



Xl4 Xl4 Xl4 

— r I/,";;-! o i ' : 

Xl4 Xl4 Xl4 

, 2k r £ 2 £ 2 2k v £ 2 s 12 

H 2k 1 -£ 2 s 12 h 2k 3 -£ 2 s 12 - 2x\^k z -£ 2 s X2 - 2xuh-£ 2 s 12 

Xl4 Xl4 

-2£ 2 Sl2 - ^ - A£ r £ 2 s l2 - 2 -h^L _ 3 ^ Si2 _ _ 2xi4£ 2 Si2 + Xi4S 2 2)£ , 

Xl4 Xl4 Xl4 

The algorithms described in the previous section are not guaranteed to yield the solutions 
in the simplest possible form; it can happen that linear combinations of solutions can be 
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factored to yield a solution of lower engineering dimension. In this case, however, the 
solutions do appear to be close to the "simplest" possible ones. 

Were we to use algorithm III, we would start by solving the equation for the irreducible 
part of E. Here, the Grobner basis has 12 vectors, giving rise to five syzygies, and three 
solutions to eq. f 14 . 3 5 1) — one with coefficients of engineering dimension two, the other two 
with coefficients of engineering dimension four. At the second stage, we have 103 basis 
tuples in the Grobner basis (limited to engineering dimension six) for E augmented by 
the auxiliary equations (I4.37P imposing reducibility. We find that that the right-hand side 
of eq. (I4.36P can be decomposed over this basis, so that ry in eq. (I4.3ip vanishes. There are 
again three solutions to the equations. At the third stage, we include the last column of E, 
corresponding to the propagator involving both loop momenta; this yields 125 tuples in the 
Grobner basis, and the right-hand side again reduces over this basis, and we end up with 
the three solutions prefigured by the solutions to the irreducible equations. In more complex 
integrals, the third stage will often impose additional constraints on the free polynomials 
obtained at the first stage, leading to fewer solutions, or more solutions with coefficients of 
higher engineering dimension. 

The two algorithms are not guaranteed to produce the same solutions, and even for the 
massless double box, they do not. The solutions will however span the same space, and yield 
the same solutions to the IBP equations for the integrals of interest. In this case, they do 
produce the same number of solutions of each engineering dimension, and the solutions are 
equivalent — the solutions produced by algorithm I can be written as linear combinations 
of those produced by algorithm III and vice versa. 

With the IBP generating vectors of eqns. (15. 8 p and (15. 9p . we can now construct IBP 
equations, 

= p £[^Ki + ^K*]> ( 5 - n ) 

where p is an irreducible polynomial in the symbols in V 22 . These equations will relate 
various nominally-irreducible P 22 to integrals with fewer propagators, but by construction 
will involve no undesired integrals. In practice, we do not need all three solutions; the first 
two suffice to produce all possible IBP equations. 

The first solution (I5.8P leads to the following equation, 

= 2Pr 2 [^i^2](A;i,A; 2 ,A ; 3,A;4)-2Pr2^4^i](A;i,A; 2 ,A;3,A ; 4) + P 2 * 1 [l](A; 1 ,A; 2 ,fc3,A;4) (5.12) 
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-2^1] (fc 3 , K h, k 2 ) + P^ilKh, k 3 , k 2 , h) - P 2 *! [l](fci, k 2 , K u ) + P£[l](k 4 , h, K 12 ) 

If we make use of the symmetries of the reduced integrals P 2l and P 2 \, or reduce these 
latter integrals in turn, this equation simplifies to, 

p;* 2 [h ■ £ 2 }{h, k 2 , k 3 , h) = P™[h ■ £i]{h, k 2 , h, h) . (5.13) 

For the all-massless double box, this equation is also a direct consequence of the symmetries 
of the integral; but the analogous statement is no longer true for double boxes with external 
masses. 

We can also use these vectors to derive equations for double boxes with more complicated 
numerator insertions, of powers or products of the basic irreducible numerators. As discussed 
earlier, we can do so by multiplying the vector by powers of invariants, which still yields a 
solution to the equations requiring that the IBP be free of doubled propagators. 

In a gauge theory, 22 irreducible double boxes can arise: 

p**[i], p*; 2 [h ■ £ 2 ], p*; 2 [(h ■ £ 2 )% p**[(h ■ £ 2 ) 3 }, p 2 **#4 • hi 
p** 2 [(h ■ £ 2 )(h ■ h)],P%[{k x ■ £ 2 ) 2 (h ■ h)}, p*; 2 [(h ■ £ 2 Wh ■ i x )i p**p 4 ■ t x )\ 
p^Kh ■ £ 2 ){k A ■ £ x )\ P**[(ki ■ £ 2 ) 2 (h ■ £,)% p*; 2 [(h ■ £ 2 f(k 4 ■ £inP^; 2 [(h ■ £,)% 
P^ih ■ £ 2 ){k A ■ £ x )% P%[(h ■ £ 2 ) 2 (k 4 ■ £i)% P^Kki ■ ^f{k A ■ £,)% P 2 **p4 • e 1 )% 

p*; 2 [(h ■ £ 2 )(h ■ £i)%P 2 M k ^ ■ ^) 2 ( k * ■ p 2%\(ki ■ ^) 4 ], p 2 *,* 2 [Oi ■ • £i)l 

P^Hh-l^ih-h) 2 }, (5.14) 

where we have omitted the momentum arguments k\, k 2 , k^, k A for brevity. In a gravitational 
theory, higher powers of the irreducible numerators may occur. 

Two of the three IBP-generating vector pairs suffice to generate all possible IBP equations 
for these integrals (the third pair yields only linear combinations of the same equations). If 
we require that the coefficients be non- vanishing in the limit e — > 0, we find 19 equations; we 
can obtain an additional equation by relaxing this constraint. This allows us to eliminate 
20 of the 22 integrals, solving for them in terms of integrals with fewer propagators and two 
irreducible master integrals, for example, 

P**[l] andP**# r 4]. (5.15) 

This reduction is the same as previously obtained with air [k^ (and presumably by oth- 
ers). Even though the coefficients may be of order e, the solutions do not involve singular 
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coefficients for the double-box master integrals. (Because we do not fully reduce the simpler 
integrals, we cannot determine whether that is also true for them.) We do not need to study 
the full set of 22 integrals; a minimal set that reduces fully with the same two IBP-generating 
vectors is, 

{P**[l] ,P**[A;r4] ,P2,* 2 Pr4) 2 ] ,P%[k<'li] ^^[h-e^-exj^^iih-e,) 2 ]} (5.16) 

In addition to the IBP equation (15. 12ft or (15 . 13f) . we have three additional equations for this 
set, arising from the first IBP-generating vector pair with prefactors k^l\ or k 2 -£2, and from 
the second IBP-generating vector pair with no prefactor. 

The form of such minimal sets will in general depend on the dimensionality of the solution 
vectors; simplifying the IBP-generating vector pairs by taking linear combinations can in 
general lead to simpler minimal sets of effectively-reducible integrals. 

The above reductions hold to all orders in the dimensional regulator e. We can also ask 
whether additional relations appear when we drop terms of 0(e) in the integrals. For the 
double box, there is only one Gram determinant which can lead to such a relation, 

G(£ 1 ,£ 2 ,1,2,A), (5.17) 

and one linear combination, 

G(£ 1 , 1, 2, A)G(£ 2 , 1, 2, 4) - G 2 J' . (5.18) 
When we use the IBP equations to reduce 

p;* 2 [G(e u e 2 , 1,2, 4)], (5.19) 

however, we find that the two irreducible master integrals (I5.15P both appear with coefficients 
of 0(e), and hence the Gram determinant fails to produce a useful relation. (More precisely, 
it provides only a relation for the divergent terms in the two integrals, but not for their finite 
terms.) The same is true for the combination of eq. (15.181) . This strongly suggests that both 
integrals (15. 15ft that are independent to all orders in e remain linearly independent when 
truncated to O(e ). (It doesn't provide a complete proof because we have not proven that 
the expressions (I5.17P and (" 15 . 1 8[) give all possible relations of this type.) 
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FIG. 6: Double boxes with external masses, with the massive legs indicated by doubled lines: (a) 
one-mass (b) short-side two-mass (c) diagonal two-mass (d) long-side two-mass. 

VI. MASSIVE DOUBLE BOXES 



In this section, we survey the IBP-generating vectors for double boxes with some of the 
external legs taken to be massive. There is one possible configuration of masses if one 
external leg is massive, as is also true if three or four external legs are massive. With 
two massive external legs, there are three possible inequivalent integrals: both massive legs 
adjacent and attached to the same loop ('short side' or 2ms); the massive legs attached to 
diagonally-opposite corners ('diagonal' or 2md); or massive legs adjacent but attached to 
different loops ('long side' or 2m/). The one- and two-mass double boxes are shown in fig. |6l 

We can use the same basis momenta and hence same form (15.1 ft and the same basic 
symbols ( 15. 2 j) as in the massless case. Following the procedure outlined in section HV] we 
find three IBP generating vectors for the one-mass double box (we take leg 1 to be the 
massive one), once again one with coefficients of engineering dimension two, and two with 
coefficients of engineering dimension four. It again suffices to use the first two vectors to 
generate all possible IBP equations; there are again 20 equations for the 22 nominally- 
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irreducible integrals, giving rise to two irreducible master integrals, say, 



P 2 **[l] andP 2 **# r l 2 ]. 



(6.1) 



The set of integrals in eq. ( I5.16P is again a minimal set that can be reduced. We will not dis- 
play the IBP-generating vectors explicitly, but they are given in a companion Mathematica 
file. 

For the long-side two-mass double box, we take legs 1 and 4 to be massive, and now 
find five IBP generating vectors, all with coefficients of engineering dimension four. There 
are again 20 equations for the 22 nominally-irreducible integrals, which we can derive using 
three of the five pairs of vectors. We can again pick the integrals in eq. (16.11) as irreducible 
masters. 

For the diagonal two-mass double box, we take legs 1 and 3 to be massive, and use 
algorithm III to find three IBP generating vectors, with the same dimensions as the massless 
and one-mass cases. Once again, we need to use only two vector pairs to generate all required 
equations, and can take the integrals in eq. ( 16. ip as irreducible masters. 

When we examine the short-side two-mass double box (taking legs 1 and 2 to be massive), 
we find our first surprise. Here we find four IBP generating vectors, all of engineering 
dimension four; but we find only 19 equations for the 22 original integrals (for which we 
need three of the four vector pairs), leaving us with three irreducible master integrals, 



The three- and four-mass double boxes lead to very complicated analytic expressions 
in intermediate stages; it is much faster (and sufficient for the unitarity approach, as dis- 
cussed earlier), to compute the IBP-generating vectors for a fixed numerical configuration 
of external momenta. We have chosen to do so; we find five IBP-generating vectors for 
the three-mass case (three with coefficients of dimension four, and two with coefficients of 
dimension six), and four vector pairs for the four-mass case (two each of dimensions four 
and six). In line with the result for the short-side two-mass double box, we find three mas- 
ter integrals for the three-mass case (which we can take to be those in eq. (16.21) ). For the 
four-mass case, we find that we need four master integrals, which we can take to be, 



p 2 ** 2 [i] , p;: 2 [h-e 2 ] , and p 2 *;#4-4] . 



(6.2) 



P 2 **[l] , Pr 2 [h-£ 2 } , P^ih-h] , and Pftih-hkfh] . 



(6.3) 



In all cases, there are no additional equations that arise from truncation to 0(e). 
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VII. THE PENTABOX 




Our next example is one of the three basic topologies that arise in five-point computations: 
the pentabox P^* 2 , shown in fig. [7J Here, we choose 2,3,5 as basis momenta, so the general 
form (14. 2 p becomes, 

v$ = cf^il + cfh^ + <^k» + cf ^ + cfhl + c®A£ , (7.1) 

where again each of the coefficients is a function of Lorentz invariants in the set of 
symbols V 32 , 

V 32 = {£ 2 1 ,£ 1 -£ 2 ,£ 2 2 ,£ 1 -k 1 ,£ 1 -k 2 ,£ 1 -k 3 ,£ 1 -k 5 ,£ 2 -k 1 ,£ 2 -k 2 ,£ 2 -k 4 ,£ 2 -k 5 ,s 12 }. (7.2) 

For this integral, we have constructed vectors both analytically and numerically; the 
numerical construction is much less memory-consuming. In both cases, the algorithms yield 
six IBP-generating vectors with coefficients of engineering dimension four, and three vectors 
with coefficients of dimension six. Their forms are too lengthy to display in the text, but are 
provided in the companion Mathematica file. There are 76 nominally-irreducible integrals 
in a gauge theory, involving powers of the three irreducible numerators, 

h-£ 2 , k 2 -£ 2 ,h-h- (7-3) 

It suffices to use the six vector pairs of dimension four to generate all possible equations for 
these integrals. We find 73 such equations, leaving us with three truly-irreducible master 
integrals, which we can choose to be, 

^; 2 [i], pfflhM ,p&[k6-ii]. (7.4) 
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Examples of these reduction equations are, 

P ^ M = 4( Xll - X23 + xa ) 

vl ^ 54 ^3*2^1 -^2] + simpler integrals , 



Pr 2 [k r £ 2 k 2 -£ 2 ] = xi 5 



+ 



Xis ~ X23 + X45 

1 + Xl5 ~ X34 ~ X45 ~ X15X45 + X23X45 + X34X45 

8(l-2e)(l-X34 -X45) 

_^ £(1 + Xl5 ~ X23 ~ X34) 

8 (1 - 2e)(Xl5 - X23 + X45)(l - X34 - X45) 

X(X15(1 - X45) + (X45 - X2 3 )(2 - X34 - 2X4 5 )))s?2^3** 2 [l] 

€(Xl5 + 2X15X23 ~ 2X23 ~ X23X34 ~ X15X45 + 2X23X45 + X34X45) 

2(l-2e)(xi5 -X23 + X4 5 ) 

2 + 2X15 ~ 3X34 ~ X15X34 + X34 ~ 2 X45 - 2X15X45 + X23X45 + 2X34X45 \ 

2(l-2e)(l-X34 -X4 5 ) / 

x Sl2 p;* 2 [k v £ 2 ] 

(1 + 2e)(l + X15 - X23 - X34)(l + X23 - X45) 



-Si2P 3 * 2 [k 



5-*l 



4(l-2e)(l-X34 -X4 5 ) 
+ simpler integrals . (7.5) 

These master integrals are independent when considered to all orders in e. Unlike the 
case of the double box, however, here we can find two linear relations between them, so 
long as we truncate at O(e ). These two relations arise from considering the following two 
integrals, 



^3,2 



g( £ ^ W k )} andP^Mi^ 1 ' 2 ' 3 ' 5 ) 
W 2 ,l,2,3,5yJ 3,2 L U 2 , 1,2, 3, 5 J 



? 2 , 1,2, 3, 57 J (7 ' 6) 
both of which are of 0(e), as discussed in section [Til El Setting the two to zero, and using the 

reductions obtained from integration by parts, we find two equations relating the masters 

in eq. (I7.4p . We can use these, for example, to eliminate the two integrals with non-trivial 

numerators in favor of P 3 * 2 [l], 

P&fa-ti] = Xl5X34X l 5Sl2 P5[l]+ simpler integrals + 0(e) , 

-Xl5 + X23 - X23X34 + X15X45 + X34X45 
p**^.^} = X15VX15(1 ~ X45) 2 + X34(l - X45)X45 ~ X23(l ~ X45 ~ X34(l + X45))) g p*» Ml 

3,2 1 2 4(1 - X34 - X45XX15 - X23 + X23X34 - X15X45 - X34X45) 12 3,2 

+ simpler integrals + 0(e) . (7.7) 

The other Gram determinants of the form suggested in section UlI El do not yield independent 
equations, but could be used instead to obtain equivalent equations. The longer denomina- 
tors in these expressions may develop poles at exceptional values of the kinematics; these 
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are presumably spurious and are softened by the behavior of the various integrals in those 
limits, but we have not checked this. (Other denominators vanish in non-adjacent collinear 
limits, for example 1 — % 34 — x 45 — > in the collinear limit k 3 || k 5 ; these are presumably 
spurious as well.) 

One may wonder whether the IBP equations are even required for reduction of the trun- 
cated integrals, given the seemingly-stronger equations arising from Gram determinants. 
However, this strength is illusory: if we use only Gram determinant equations (including the 
IBP-like ones built from Gram determinants of the form given in eq. (I3.23P ). we find only 24 
equations for 29 of the 35 integrals with numerators of dimension eight or less. (The IBP-like 
determinants provide two of these equations.) This would leave five seemingly-irreducible 
integrals as masters; of course, using the IBP equations, we could then reduce all of these 
to the scalar integral P 32 [l]. If we consider the complete set of 76 integrals, we see another 
problem with using Gram-determinant equations alone: 20 of the integrals (those with four 
powers of £2) are ultraviolet divergent, which prevents us from using these equations to sim- 
plify them. In addition, even amongst the ultraviolet-finite integrals, we are left with five 
master integrals (there are 51 equations in total that we could derive). 

VIII. A SIX-POINT EXAMPLE 



3 4 




1 6 

FIG. 8: The double pentagon P 3 ** 3 . 

We will consider one example of a six-point integral, the so-called double pentagon P 33 , 
shown in fig. [HJ In contrast to the pentabox P 3 *2 considered in the previous section, we 
find that this integral can be reduced to simpler integrals entirely using Gram determinant 
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equations alone. Indeed, not only can integrals with non-trivial numerators be reduced, 
but the scalar integral itself, Pg 3 [l], can also be expressed in terms of simpler integrals 
(pentaboxes and products of one-loop pentagons, themselves reducible) via algebraic identi- 
ties. We don't even need those Gram determinants equivalent to IBP equations to perform 
these reductions. 

The double pentagon has two irreducible numerators, which we can pick to be k§ ■ £\ and 
k\ ■ £ 2 . There are thus 33 formally-irreducible integrals that arise in a gauge theory. The 
first thing to notice is that they are all ultraviolet-finite, so one of the obstructions that 
existed in the pentabox case to use of Gram-determinant equations alone for a complete 
reduction is absent here. There are 15 integrals of engineering dimension eight or less; 
exclude -P3*3[(^6'4) 4 ] an d -R^K^r^) 4 ], and examine the remaining 13 integrals. We can find 
13 independent equations for them by starting with the following identities, 



(9(e) 



r 3,3 



P G 



?i,l,2,3,6 
'2,1,2,3,6 



with prefactors p = 1, ki-£ 2 , (k 1 -£ 2 ) 2 , k 6 -£i, (£; 6 -4) 2 , k\-£ 2 k { 



5.1) 



(9(e) 



-"3,3 



P G 



?i,4,l,2,3 
?i,4,4,5,6 



with prefactors p = l,ki-£ 2 , (k 1 -£ 2 ) 2 , k 6 -£i, (k e -£i) 2 ; and 



(8.2) 



(9(e) 



r 3,3 



G 



? i,^2,l,2,4 
?i,4,3,5,6 



and (9(e) 



-"3,3 



G 



4,4,1,2,5 
4,4,3,4,6 



(8.3) 



Similarly, if we examine the full set of 33 formally-irreducible integrals, we find 33 inde- 
pendent equations. 



IX. CONNECTION TO GENERALIZED UNITARITY 

In this section, we use generalized unitarity to give a heuristic explanation for the struc- 
ture of the results presented in previous sections. We begin, as in section [Til with a discussion 
at one loop. 

In basic unitarity at one loop, we examine the branch cut of amplitudes channel by chan- 
nel. In each channel, the branch is a phase-space integral over a product of tree amplitudes. 
In the present paper, we are focused on loop integrals rather than complete amplitudes, so 



the equivalent statement — dating back to the Cutkosky rules 
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52| of the 1960s — is the 



expression of the branch cut in terms of phase space integrals of scalar tree diagrams. The 
ordinary cut may be obtained by cutting two propagators, that is replacing the propagators 
by positive-energy delta functions which put the intermediate state on shell, 

► 2tt(5 (+) ((£ - Kf) . (9.1) 



(£ - Kf + is 

(The (+) superscript indicates the restriction to positive energies.) There is nothing stopping 
us, however, from cutting more than two propagators, and this is the idea behind generalized 
unitarity. The solutions to the delta function constraints will in general then be complex, and 
so the delta functions must be understood in a more general sense, as contour integrals with 
the contours chosen to encircle the common solutions to the constraint equations. The idea 
of generalized unitarity was first applied as a practical tool for computation of amplitudes 



by Bern, Dixon, and one of the authors 



23j ]. It was later combined with the use of complex 



momenta to give a general algebraic solution to finding the coefficients of box integrals 171 ]. 
and used to derive a g eneral and numerically-applicable technique for triangle and bubble 



integrals by Forde 



19|. 



If one cuts as many propagators as possible, one arrives at maximal unitarity, as used for 



example in ref. 
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53]. In old-fashioned language, this is equivalent to looking for 'leading 



singularities', discussed in a modern incarnation in ref. 54]. 

But how many propagators can we cut? If we examine a one-loop amplitude with all 
external momenta taken to be massive (so that infrared singularities are tamed), we can 
take the dimensional regulator e to zero, and perform the integrals in four dimensions. 
(Ignore the ultraviolet-divergent bubble in this discussion.) Each delta function will impose 
one constraint; because we have four components, we can have up to four delta functions. 
Attempting to impose additional delta functions will in general yield no solutions. (More 
precisely, because we will have more delta functions than integrals, the result will itself be 
a delta function rather than an ordinary function.) This in turn implies that functions with 
additional propagators may be determined in terms of functions with up to four propagators, 
as there are no additional degrees of freedom. In this case, all pentagons or higher-point 
integrals are reducible to sums of boxes and lower-point integrals. 

The generalization of this observation to higher loops is straightforward. At each loop 
order, we have an additional four components. We can thus cut an additional four propaga- 
tors. When considering infrared- and ultraviolet-finite integrals, then, we expect that only 
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those with up to four propagators per loop momentum will be algebraically independent. 
At two loops, this means that integrals with more than eight propagators, or more than four 
propagators involving a single loop momentum, will be reducible into simpler integrals. 

Of course, the integrals of interest are in general infrared divergent. While the loop 
momentum is formally D-dimensional, so long as we keep the external momenta in four 
dimensions, the additional components \i can only appear at one loop as /i 2 , on which we 
can impose one additional delta function. Thus when considering one-loop integrals to all 
orders in e, the pentagon integral must be taken as an additional independent integral, while 
higher-point integrals remain reducible. Now, the algebraic independence of the pentagon 
only manifests itself at 0(e); terms through O(e ) are reducible to sums of boxes. This 
reducibility is not manifest in our heuristic discussion; but it suggests the conjecture that 
the reducibility of integrals to O(e ) follows the pattern of massive reductions. 

What happens at two loops? We now have two loop momenta, and correspondingly two 
e-dimensional vectors, fi\ and /12. These can now appear in integrals in the form of three 
independent quantities, fif, and /ii ■ \i2- We could impose additional delta functions 
on each, corresponding to cutting three additional propagators. For two-loop diagrams, 
we therefore expect any integral with more than eleven propagators, or more than five 
propagators involving a lone loop momentum, to be reducible to all orders in e. This is 
exactly what we found in section IIII1 

Different propagators lead to different branch points (or branch surfaces, for the many- 
complex- variables functions we are considering). Accordingly, the algebraic independence 
of uncut propagators is clear. The algebraic independence of non-trivial numerators is less 
clear, as one might imagine algebraic relations between them. (Indeed, there are clearly 
algebraic relations between different powers of numerators, as seen in reduction equations 
elsewhere in the literature or in previous sections.) Heuristically, we do at least expect an 
upper bound, 

# irreducible non-trivial numerator integrals < 11 — ^propagators . (9-2) 

This bound is respected by the explicit results in previous sections. 

The question of algebraic independence when truncating integrals to 0(e) is more sub- 
tle. If we adopt the conjecture above suggested by one-loop results, it would imply that 
only truncated integrals with up to eight propagators (and up to four involving each loop 
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momentum) are algebraically independent. This is in agreement with the reducibility of the 
double-pentagon P 3 *g discussed in section IVIIII We might be further tempted to conjecture 
that the number of independent integrals with non-trivial numerators is limited to eight less 
the number of propagators. This would imply that there are no independent pentaboxes 
with irreducible numerators, which is in fact true. Thus the bound is respected by the 
pentabox results discussed in section IVHt and also by the results for some double boxes; but 
it is violated for the short-side two-mass double box, as well as for double boxes with three 
or four external masses. The precise manner in which the heuristic picture breaks down 
remains to be clarified. 

X. CONCLUSIONS 

Knowledge of an integral basis plays an important role in modern unitarity calculations. 
In this paper, we have given an outline of a basis for planar integrals (with massless prop- 
agators) at two loops. We distinguish two kinds of bases: the first, a set of integrals which 
are linearly independent to all orders in the dimensional regulator e; the second, in which 
linear independence is required only through O(e ). 

Smirnov and Petukhov [29{] have recently shown that the integral basis resulting from 
integration by parts is finite. We have delineated an explicit finite set of integrals which 
contains a minimal basis, and given an explicit procedure reducing an arbitrary planar two- 
loop gauge-theory integral to an element of this set. The set contains only integrals with four 
or fewer external legs attached to each line in the vacuum graph. All irreducible numerators, 
whose number depends on the external legs, are allowed in this finite set of integrals. The 
final basis contain only a subset of these integrals. 

In order to reduce the set further, we then introduced an approach to generating 
integration-by-parts equations which involve only integrals in the desired set, along with 
simpler integrals (with propagators omitted), and avoiding integrals which are not ordinary 
Feynman integrals. For each of the integrals in the above set, one can solve for these vec- 
tors, and then determine the set of independent master integrals that make up the first, 
D-dimensional, basis. Unlike the situation at one loop, the reductions, and more impor- 
tantly, the number of independent integrals, depend on the masses of the external legs. We 
gave a few examples of this type of calculation, but leave a complete study of the integrals 
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to future work. 

We also introduced a special set of numerators, built using Gram determinants, which 
provide equations that yield identities for integrals truncated to O(e ). These equations 
reduce the D-dimensional set of independent master integrals to a smaller set making up 
the second, "regulated four-dimensional" basis. 

The general arguments as well as the notion of IBP-generating vectors and additional 
0(e) identities generalize to non-planar integrals as well as to higher- loop integrals. We also 
expect them to generalize from the massless propagators considered here to integrals with 
massive propagators. It would also be interesting to explore the analog of the (9(e) identities 
for integrals in two- and three-dimensional field theories. We gave a heuristic argument for 
understanding the basis in terms of generalized unitarity; it would be interesting if it could 
be developed further to a complementary derivation for the reduction of integrals with 
formally-irreducible numerators to an independent set of master integrals. 

The defining equations (14.91) for the IBP-generating vectors can also be thought of as 
defining a 'surface' or variety in the space whose coordinates are given by the different 
monomials in V ( 14.4ft . It would be interesting to explore its connection with the Grass- 
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551 ] introduced in recent explorations of integral coefficients in the M = 4 theory. 



The integral basis appropriate for the M = 4 theory should presumably make manifest (up 
raxed divergences) its extended symmetries (conformal and dual conformal symme- 
56j|). and may make natural use of twistorial integrands such as those discussed in 



to in 
tries 
refs. 



37, 
57]. 



Our approach to solving the required equations ( 14. ip for the IBP-generating made use of 
Grobner bases, and in particular the standard Buchberger algorithm 43|, |44| for computing 
them. The present implementation of the algorithm (coded in Mathematica) performs well 
for simple cases like the double box, but slows down and requires large memory in its 
intermediate stages for integrals with more legs or many massive legs. It would be worthwhile 
exploring the use of more modern algorithms, such as those of Faugere [50;] for computing 
the required Grobner bases. 
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